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I. INTRODUCTION , 

Wide-angle celestial structures, such as bright comet tails and nearby 
galaxies and clusters of galaxies, rely on photographic methods for quantified 
morphology and photometry, primarily because electronic devices with comparable 
resolution and sky coverage are beyond current technological capability. CCDs with 
pixel numbers measured in the millions have become commonplace, but fail by orders 
of magnitude in providing sky coverage for large structures with good resolution. 
Even computers are daunted by the task of processing bits of information numbered in 
the thousands of millions. The photographic process succeeds for all Its obdurate 
properties, and a picture speaks a thousand megawords. 

Nevertheless, seeing is not quite believing, because photographic density is 
not simply related to luminous intensity. Morphology is clearly dependent on 
photometry, particularly in non-linear parts of the response, and astrophysically 
important numbers are dependent on both,__ While stars of known magnitude can be used 
to calibrate images of other stars, their use in determining the surface brightness- 
es of extended images is not straightforward. 

Brandt (1985) has emphasized that stars of known magnitude are natural 
calibrators that should lead to accurate calibration of photographic plates if a 
suitable theory is forthcoming, because In principle all the necessary information 
resides In the exposed emulsion. Rather than being an unavoidable hindrance to 
isophotometry, star images would play a pivotal role in calibration. While spot 
sensitometry is the canonical means of plate calibration, it is not always available 
on plates of interest, and suffers from its own intrinsic calibration problems. 

The present work is an examination, more-or-less from first principles, of 
the problem of the photometry of extended structures and of how this problem may be 
overcome through calibration by photometric standard stars. Previous work is 
summarized In 5TII. In IfHI the perfect properties of the ideal f ield-of-view are 
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stated in the guise of a radiometric paraxial approximation, in the fond hope that 
f ields— of —view of actual telescopes will conform. Fundamental radiometric concepts 
are worked through in 1TXV before the issue of atmospheric attenuation is addressed 
in W. The independence of observed atmospheric extinction and surface brightness 
leads off the quest for formal solutions to the problem of surface photometry in 
1TVI, 

Methods and problems of solution are discussed, in WII. The development 
relies on the powerful method developed by Zou, Chen, and Peterson (1981), but puts 
it on a radiometrically sound footing by solving the problem of the radiometric 
zero-point, and by ensuring through conservation of energy that atmospheric extinc- 
tion and the radiometrically related sky brightness are actually part of the 
solution. Standard star image profiles need not conform to any theoretical ideal, 
but images must not be too distorted (see §7*6); on all counts, the present 
formalism is comparatively assumption— free. In its final form, all quantities that 
appear are dimensionless as they should be. 

liller (1985) has raised the spectre of reciprocity failure (RF), and 1TVII 
does not shrink from it* Rather, the spectre is confronted in the spirit of 
standard stars and shown to be chimerical in that light, provided certain rituals 
are adopted. RF (and a good deal else besides) can be exorcised in future studies 
of comets provided that standard stars are inserted just as sensitometry spots are 
inserted . 

After a brief discussion in Will of Baker-Sampson polynomials and the 
vexing issue of saturation, 1TIX embarks on a pursuit of actual numbers to be 
expected in real cases. While the numbers crunched are gathered ex nihilo, they 
demonstrate the feasibilty of Newton’s method in the solution of this overdetermin- 
ed, nonlinear, least square, multiparametric , photometric problem. This is no small 
feat, and should encourage future research. 
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1TX concludes the Report with a summary of Reports I and II (Dorband et al. 
1985; Dorband and Usher 1985); the definition of the search area for photometric 
standards in the P/Halley path, and their incidence and usefulness in future 
photometric reduction, is described. 

Errors of one sort or another are bound to exist in spite of my best 
intentions, and I apologize for them in advance; I would appreciate your communicat- 


ing them to me. 
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II. OTHER METHODS . 

The first step in the process of converting density to surface brightness is 
the digitization of the image (see e.g. Klinglesmith 1983). Spot size should be 
larger than the scattering radius and grain separation, but large enough to reduce 
the noise of the output. Stock and Williams (1962) review 3 methods of data 
acquisition: (i) schraf f ierkassette , (ii) Fabry photometry, and (iii) extrafocal 
images, and they discuss the advantages and shortcomings of each. In (i), a moving 
image introduces the Intermittency Effect (combined with RF too if images are both 
point and extended). In (ii), the image size of point and extended images is the 
same, but complex lens systems are needed. In (iii), extrafocal images destroy 
morphological information, but Stock and Williams make an interesting point: 
extrafocal images are useful calibration tools when imposed in a separate exposure 
of the same length on a portion of the plate different from an extended object of 
Interest. This concept of separate standard star sequences Is a precursor to our 
solution to the problem of RF for trailed standards on comet plates (§7.6), except 
that. In light of recent advances (e.g. ZCP, Klinglesmith and Rupp 1984; Warnock and 
Klinglesmith 1984, and this Report), there is no need for the images to be out of 
focus! Moreover, there is the added bonus that atmospheric extinction can be found 

as part of the solution (§6.3 et seq.). 

Stock and Williams (1962) also discuss sensitometric spot and wedge calibra- 
tion. This is far and away the most common method nowadays, but it has drawbacks 
(see Agnelli et al. 1979); the object and the calibrators are not simultaneously 
exposed, and ambient conditions are oftentimes different; also, the characteristic 
curve may differ from one part of the plate to another. In our proposed method, 
ambient conditions will be different only if outside conditions are different, and 
these are amenable to experimental control and are capable of being monitored. As 
for not having the H&D (or characteristic curve) the same from one part of the plate 
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to another, multiple exposures are necessary and desirable in any case. In fact, 
standard stars can be used to calibrate the wedges. 

The methods of Kormendy (1973), Agnelli et al, (1979), and Feitzinger et al. 
(1983) depend on assumptions regarding the profiles of field stars. Image symmetry 
is desirable and assumed to be independent of apparent magnitude, Kormendy (1973) 
determines the magnitude zero-point by simultaneous photoelectric photometry of the 
sky. The methods are limited by their assumptions and complexity, and exiguous in 
their lack of simultaneous solutions for sky brightness and extinction. Owing to 
the crucial role played by the sky brightness in Type II detections and thus in the 
definition of morphological structure (e.g. Sandage 1972), it Is essential that the 
faint end of the H&D curve be well defined, but Feitzinger et al. (1983) show that 
this Is not easily accomplished, 

Zou, Chen, & Peterson (1981) [herein ZCP] give a method for star profiles of 
arbitrary shape, and so allow somewhat for a less-than-perfect world. In the 
exemplification of their method, the zero-point of the calibration (which Is de 
facto the background sky brightness) is assumed, and there is no allowance for 
atmospheric extinction. The methods of ZCP and Agnelli et al. (1979) have been 
implemented by Klinglesmith & Rupp (1984) and Warnock & Klinglesmith (1984), though 
the problem of the photometric zero-point, extinction, and sky brightness, remain. 
The great promise shown by the moment-sum method of ZCP has suggested to the 
Steering Committee of the IH-W LSPN that its applicability to P/Halley, and future 
comet apparitions, be examined. 
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III. THE RADIOMETRIC PARAXIAL APPROXIMATION. 

Wide-field imagery is accomplished normally with telescopes of the classic 
Schmidt design. Such cameras are stopped at the center of curvature of a concave 
spherical mirror, and so produce images that are free of coma and astigmatism. But 
they are afflicted with spherical aberration, and this is reduced to fifth order by 
an aspheric corrector situated at the camera stop. Thus each image lies on an 
optical axis that passes through the camera stop, and on a spherical focal surface 
whose center is at the center of curvature of the primary mirror. The basic 
principle of the Schmidt system therefore is that the telescope has no preferred 
optical axis, and so yields excellent Images to high order in all field directions 
(Bowen 1960). 

Distance-to-center effects will arise however because the corrector has a 
preferred axis. Images far from the axis of optical symmetry will suffer losses In 
Irradiance owing to vignetting. There is also a deterioriation of Image quality 
owing to residual aberrations; for oblique Incidence, the residual fifth-order 
spherical aberration sets the limit on the f ield-of-view for good images, and coma 
and astigmatism also develop owing to the relative tilt of the corrector. 

In general, the Radiometric Paraxial Approximation (the RPA) is valid for 
those directions sufficiently close to the axis of optical symmetry that losses In 
irradiance are negligible. For example, the unvignetting field of the 48-Inch 
Palomar Schmidt camera Is about 5.4 degrees; even at the corners of the 14-inch 
square photographic plates, the loss of irradiance is estimated to be less than 0.2 
mag. (Minkowski and Abell 1963). 

In the RPA, aberrations affect surface photometry only through a loss of 
morphological discrimination, since radiometric theory accounts for all the energy. 
The photometry of point sources is similarly unaffected. A radiometric calibration 
of the entire FOV of the telescope can be accomplished by standard stars sufficient 
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ly close to the direction of interest. An accurate morphological-photometric 
description of structures larger than the field either of negligible aberrations or 
negligible vignetting may in theory be achieved with the help of the spatial 
response of the camera to point sources, but this is beyond the scope of the present 
treatment. In addition, chromatic aberration is a fundamental limitation to the 
morphological discrimination of any camera with a refracting element, and is present 


even in the RPA, 
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IV. RADIOMETRY 
§4.1 PIXEL RADIANCE 

Let L(A ) be the Spectral Radiance of a source at infinity incident upon 
a telescope of nominal aperture diameter D, focal length f, and effective 
collecting area A g - VirGD 2 . The function G accounts for geometrical effects such 
as vignetting and shadowing of the primary. 

A pixel area a - f 2 fi in the focal plane subtends a solid angle ft of the 
sky. The radiance L(A) suffers losses owing to the transmittances T t (A) and 
T^(A) of the telescope and filter respectively. If the detector has a quantum 
efficiency R^(A), then the System Response is: 

(4.1.1) S(A ) - T t (A)*T f (A)-R d (A). 

The detected spectral radiance is: 

(4.1.2) F(A) - L(A ) S(A), 

and the detected Spectral Illumination is: 

(4.1.3) 1(A) = / F(A) dw. 

Q 

The heterochromatic Illumination is found by integration over the bandpass of 
the system response; this bandpass is symbolically represented by S* Thus: 

(4.1.4) I(S) * / 1(A) dA - / / F(A) dA d ai . 

e b n 

The heterochromatic Power detected is: 

(4.1.5) P(e) - A 1(e) -A J J F(A) dA d<y, 

esi 

assuming A g to be independent of A and on The total Energy detected over 
integration time T is: 

Q(e) = A / 1(e) dt * A { / / F(A) dt dA d w . 

T Ten 


(4.1.6) 
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Each integration in Equation (4.1.6) entails loss of information; in 
self-explanatory notation, appropriate mean quantities are: 

(4.1.7) <H> - J H(x)dx / J dx 

X X 

where H is some function and x might be solid angle, time, or wavelength. Thus 
increased temporal, spatial, and spectral resolution, will result from 
decreased Integration time, pixel size, and bandpass. 

The heterochromatic brightness in magnitudes of solid angle & is: 

(4.1.8) m(B) = m (o) (8) - 2.5 log 1(8), 

where m^°^(B) is the magnitude zero-point appropriate to bandpass 8. Mono- and 
heterochromatic quantities are: 


Spectral Radiances 

L(X),F(X) 

erg/cm 2 / sec/ s ter /A, 

Spectral Illumination 

I(X) 

erg/cm 2 /sec/A; 

Illumination 

KB) 

erg/cm 2 /sec. 

Power 

P(6) 

erg/sec. 

Energy 

Q(6) 

erg. 

Brightness 

m(B) 

mag. 
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§4.2 HETEROCHROMATIC SYSTEM PARAMETERS 
Broadband filters deliver heterochromatic measures of the radiant 
intensity. Thus it is necessary to characterize the incident spectral radiance 
L(X) and system response S(X) with the help of suitably defined heterochromatic 
parameters. These may be divided into three categories: 

Case (i): The "System Response Solution" gives approximate brightnesses at a 

precise wavelength X^: 

X^ is the Mean Wavelength, precisely defined solely by the System Response 
S(X). Thus in this case, a precise knowledge of the characteristic wavelength 
X^ is offset by uncertainty in brightness. 

Case (ii): The "Isophotal Solution" gives precise isophotes at an approximate 

wavelength X^: 

X is the Isophotal Wavelength, defined by the condition that energy be 
precisely conserved. Thus isophotes are precise, despite uncertainty In the 
wavelength at which energy is conserved. 

Case (iii): The "Effective Wavelength Solution" gives approximate brightnesses 

at an approximate wavelength X^: 

X^ is a weighted mean as in case (I), except that in this case the weighting 

function is the entire radiant throughput S(X)L(X). 

Unlike X , both X and X make concessions to the incident spectral 
o’ I e 

radiance. In all cases, uncertainties are necessitated by the broadband nature 
of the problem. In addition, heterochromatic measures may be weighted either 
by photon energy (§4.3) or pulse counts (§4.4). 
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§4.3 ENERGY INTEGRATION 
Following equation (4.1.4), let the function: 

(4.3.1) F(A) = S(A)L(A) 

be a measure of the spectral radiance delivered by the detector. Since L(A) 
has units of energy, equation (4.3.1) is the appropriate function for energy- 
sensitive devices when integrated over wavelength. The heterochromatic rad- 
iance F Is thus: 

(4.3.2) F - J F(X)dX - / S(A)L(A)dA . 

S(A) is determinable in the laboratory. F is measurable, being the observed 
radiance weighted by the system response. The functions L(A) and F(A) are 
generally unknowable through broadband photometry, but It is still possible to 
gain some knowledge of them by approximate means. 

Since L(A) is unknown, it may be approximated by a Taylor expansion 
about some fiducial wavelength A: 

(4.3.3) L( A) = L(A) + L f (A) (A-X) + L"(A) (A-X) a /2! + ... 

where It is assumed that the spectral gradients L f , L”, . ..» are well-behaved. 
The fiducial wavelength A in equation (4.3.3) is free to be chosen so as to 
facilitate the solution. 

One way of generating Cases (i-iii) is as follows: 

(i): The more terms retained in equation (4.3.3), the more accurate the 
representaion of L(A); but L(A) and its gradients are generally unknown. 
Nevertheless, the disposable parameter A can be chosen to ensure that at least 
one term of equation (4.3.3) be precisely accounted for; clearly this term 
should be the one that nullifies the greatest source of uncertainty in the 
series expansion. Thus if the contribution of the second term in equation 



§4.3 p 2 

(4.3.3) were to be nullified by a suitable definition of then L(A) will be 

accurate automatically to 2nd order. The resulting characteristic wavelength 

A will be an average value of A weighted by the System Response Function. A 
o ° 

is the Mean Wavelength. 

(ii) : By the Mean Value Theorem, some wavelength A exists such that: 

(4.3.4) F - L(X ± ) / S(A) dA 

precisely. In this case, only the first term in equation (4.3.3) is retained, 
but the wavelength A = A must be sought by approximation through the series 
expansion. A^ is the Isophotal Wavelength. 

(iii) : Alternatively, X might be sought by the condition that the Effective 
Wavelength A^ be a weighted mean that accounts for both the incident spectral 
radiance and the losses of the optical path. 

These parameters are derived next. 

Case (i). 

It follows from equations (4.3.2) and (4.3.3) that: 

(4.3.5) F - L(*) / S(A)dA + L’(*) / S(A)(A-X)dA 

+ %L"(*) / S(A)(A-X) 2 +... 

Since knowledge of L', L", ... is generally imprecise or lacking, L(^) will be 
most accuratly determined through the measured quantity F if we specify A 
through the condition: 

(4.3.6) / S(A)(A-X)dA - 0. 

The quantity A^ so specified is the Mean Wavelength: 

(4.3.7) A q - / S(A)AdA / / S(A)dA . 

Since A q is weighted by the Response Function, it lies somewhere in the 

middle of the range of transmitted radiances. If S(A) is symmetric in A, then 
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A coincides with the wavelength of maximum response S (A ). The 3rd term of 
o mm 

equation (4.3.5) is a measure of the size of this bandpass; it yields the 
second order moment: 

(4.3.8) y 3 (A Q ) - J S(A)(A-A o )*dA / / S(A)dA . 

Higher order moments define the shape of the bandpass. 

Substitution of equations (4.3,6) and (4,3.7) in (4.3.5) gives: 

(4.3.9) F - [L(A q ) + ^P 2 L"( A q ) + ...] / S(A)dA . 

Equations (4.3.7), (4.3.8), and (4.3.9), comprise the System Response Solution 
to the heterochromatic problem. 

Case (ii). 

By equation (4.3.4), X^ must be chosen to ensure that: 

(4.3.10) F - L(A i ) / S(A)dA . 

But equation (4.3.9) is also true; thus: 

(4.3.11) L(A ± ) - L(A q ) + ^P 2 L"(A q ) + ... 

Moreover, equation (4.3.3) is also true for any fiducial wavelength, say A . 
Thus: 

(4.3.12) L(A) - L(A q ) + l'(A o )(A-A q ) + ^L"( X 0 )( A - X 0 ) 2 + . . . 

Setting X - X^ and comparing equations (4,3.11) and (4.3.12), it follows that 
the Isophotal Wavelength is: 

(4.3.13) A - A + ^P 2 L M (A )/L'(A ) + ... 

i o o o 

Equations (4.3.10) and (4.3.13) comprise the Isophotal Wavelength Solution to 
the heterochromatic problem. 

Case (iii). 

Instead of equation (4.3.7), the weighting function may be changed to 
give the Effective Wavelength: 
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(4.3.14a) X g = / AL(A)S( A)dA / / L(A)S(A)dA 

(4.3.14b) ■ A q + y 2 L'(A o ) / L(A q ) + ... 

with the help of equations (4.3.7) and (4.3.12). Equations (4.3.9) and 
(4.3.14) comprise the Effective Wavelength Solution to the heterochromatic 
problem. The cases are summarized in Table 4-1. 
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§4.4 PULSE COUNTING 
The photon spectral radiance: 

(4.4,1) N(A) =* AL(A)/hc 

gives rise to the photon radiant intensity: 

(4.4.2a) H - A e J N(A)S(A)dA 

(4.4.2b) - (A g /hc) / AL(A)S(A)dA, 

where h is Planck's constant and c the speed of light. The function: 

(4.4.3) f - J N(A)S(A)dA 

is the pulse counting analogue of equation (4.3.2). As in equation (4.3.3), 
£(X) may be expanded in a Taylor series about the fiducial wavelength X to 
give an equation analogous to (4.3.5): 

(4.4.4) f - N(X) J S(A)dA + N'(X) / S(A)(A-£)dA 

+ ^N" (X) J S(A)(A-X) a dA +... 

The arguments of §4.3 for L(A) also apply to N(A). The results are listed in 
Table 4-2, where the relation (4.4.1) gives: 

(4.4.5) he N(A q ) - A q L(A q ); 

(4.4.6) he N'(A o ) - A q L'(A o ) + L(A q ); 

(4.4.7) he N” (A ) - A L" ( A ) + 2 L’(A ) . 

' o o o o 
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§4.5 FWHM BANDWIDTH 

The integrated System Response JS(X)dX is fundamental to all hetero- 
chromatic expressions such as occur in equations (4.3.7) to (4.3.10) and in 
Tables 4-1 and 4-2. For the sake of simplicity, define S by: 

(4.5.1) 8S = Js(X)dX , 

where : 

(4.5.2) 6 = A 2 " X 1 

Is the bandwidth equal to the full width at half maximum (FWHM) of S(X), 
regardless of weighting. Thus X^ and are defined by: 

(4.5.3) S(X X ) = S(X 2 ) = ^(XJ 

where X is the wavelength at which S(X) reaches the maximum S^. Equation 
(4.5.1) allows the integrated system response JS(X)dX to be replaced symbol- 
ically by the simpler expression 8S. 

An approximate formula for 8 may be derived from the fact that the 
gradient of S(X) is zero at X^. Thus, correct to second order: 

(4.5.4) S(X) “ S m + ^(X-XJ 2 + ... 

Equations (4.5.3) and (4.5.4) give: 

(4.5.5) X. - X ± (-S /S")* + ... (k=l ,2) , 

k m m m 

so that: 

= 2(-S /S”)^ + ... 
mm 


(4.5.6) 


8 
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§4.6 

SPECIAL SYSTEM RESPONSES 

Square, quadratic, or 

cubic, responses S(X) (a<X<b) give the following 

results: 



Square: 



(4.6.1a) 

H 

-'C 
s— ' 
to 

K, (a<X<b) , 

(4.6.1b) 

m 

0, (a>X>b) . 

(4.6.2) 

X 

0 

( b+a ) / 2 , 

(4.6.3) 

M 2 - 

(b-a) 2 /12. 

(4.6.4) 

6 - 

(b-a). 

(4.6.5) 

S - 

K. 

Quadratic : 



(4.6.6) 

S(X) = 

K( X-a) (b-X) , (a<X<b) . 

(4.6.7) 

X = 
0 

(b+a)/2, 

(4.6.8) 

U 2 « 

(b-a) 2 /20, 

(4.6.9) 

0 - 

(b-a)//2. 

(4.6.10) 

s - 

K(b-a) a /3/2. 

S(X) is symmetric and 

the mean wavelength is halfway between the short and long 

wavelength cutoffs. 



Cubic : 



(4.6.11a) 

S(X) “ 

K(X-a)(b-X)(X- c ) , (a<X<b), (c<a). 

(4.6.11b) 

- 

K(X-a)(X-b) (X-c) , (a<X<b) , (c>a). 

(4.6.12) 

X - 

O 

(b+a)/2 - (b-a) 2 /20c + .... (|c|»b). 

(4.6.13) 

y 2 - 

[1 - (b+a) /2c + ...] (b-a) 2 /20, (|c|»b). 

Exact expressions obtain for 

all functions in the cubic case, but first order 


expressions in parameter c show up the effect of skewness in the response. The 


§4.6 p 2 


mean wavelength is displaced from the halfway mark by an amount that depends on 
the sign and magnitude of c. In the event that |c| * ", the cubic reduces to 
the quadratic case with appropriate re-definition of K. 

Results for asymmetric Poisson or the symmetric Gaussian distribution 
cannot be expressed in closed form owing to finite (if sometimes small) 
probabilities of a positive response at ’’negative” wavelengths; nevertheless 
the Gaussian case is addressed by Golay (1973), who gives an approximate half- 
width of 2.36y about the mean wavelength. This compares to the the exact 
Gaussian half-width (Bevington 1969) of Y « 2o /21n2 = 2.355(7, where a is the 

standard deviation of the distribution. 
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§4.7 MEAN PIXEL ENERGIES 
Equations (4.1.6) and (4.3.1) give: 

(4.7.1) Q, - A / / / F(A) dt do< dA 

J e T 0) A 

J 

on the assumption that A g is a constant. The measured energy Q is related to 
the incident radiance through equation (4.3.2), viz. 

(4.7.2) /f( A)dA = Js(X)L(A)dA - F 

where F may be chosen from Table 4-1 and equations (4.3.9) and (4.3.10). Thus 
in cases (i) and (iii): 

(4.7.3) F - [L( A q ) + ^y 2 L"(A q ) +...] / S(A)d*, 
while in case (ii): 

(4.7.4) F = L(*.) / S(A)dA. 

where the arguments of L are given by equations (4.3.7) or (4.3.13). 

In all cases, the ubiquitous integrated system response may be express- 
ed in terms of S and the FWHM bandwidth 3 according to equation (4.5.1). So 
equations (4.7.3) and (4.7.4) are more simply written either as: 

(4.7.5) F = 3 S L(* o ), 

[where the argument of L identifies the relevant case (i) or (iii)], or as: 

(4.7.6) F = 3 S L(A i ) 
in case (ii) . 

On the assumption that 3 and S are independent of time and direction, 
equation (4.7.1) becomes: 

(4.7.7) Q » A 3 S / / L(X) dt d&> 

J e T Cd 

j 

where the characteristic wavelength ^ is summoned to represent either of the 
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arguments of L in equations (4.7.5) and (4.7.6). 

From Table 4-1, it is clear that L(X) may depend on one or both of 
time and direction through one or both of L and J; i.e. the spectral radiance 
L(X) may change in intensity or shape over time or direction. Mean values are 
therefore called for, in the sense of equation (4.1.7). Thus equation (4.7.7) 
may be written as: 


(4.7.8) Qj - A g T 8 S u. 

where it goes without saying that the abbreviated notation Lj represents an 
average both over time and the solid angle of pixel j, and is to be evaluated 
according to the prescriptions and wavelengths characteristic of cases (i), 
(ii), or (iii) [Table 4-1]. 

The radiance at the aperture of the telescope may arise from the 
emissivity both of a source "s" as well as an intervening medium "m" . Also, 
the detector "d" itself will contribute to the detected energy. Thus equation 
(4.7.8) leads to the measured total energy for pixel j: 

(4.7.9a) Qj(s+m+d) = Q ^ C s ) + Qj(m) + Qj(d) 


A c T e S ^ [L j (s) + L j (m) ] + Q j (d) - 


(4.7.9b) 

The radiances L(s) and L(m) are effectively unrelated and are thus linearly 
additive, and each has its own characteristic wavelength. 



Table 4-1. 


Optical System Parameters for Energy Integration. 

BRIGHTNESSES: 

Cases (i), (iii) 

F =■ [L(A q ) + h\i a L”(A ) +...] J S(A)dA 
WAVELENGTHS : 

Case (i) 

A q = / S(A)AdA / J S(A)dA 
V a m J S(A)(A-A Q ) a dA / / S(A)dA 
Case (iii) 

A - A + v 2 L'(A )/L(A ) + ... 
e o oo 

Table 4-2. 

Optical System Parameters for Photon Counting. 

BRIGHTNESSES: 

Cases (i), (iii) Case (ii) 

f - [N(A q ) + ^y 2 N"(A ) + ...] / S(A)dA f - N(A 1 ) J S(A)dA 

WAVELENGTHS: 

Case (i) Case (ii) 

A q « J S(A)AdA / / S(A)dA X i “ A 0 + ^ y2N "( X 0 ) 1 N ’( A 0 ) + ••• 

y 2 = / S(A)(A-A Q ) 2 dA / / S(A)dA 
Case (iii) 

A = A + y 2 N'(A )/N(A ) + ... 
e o o o 


Case (ii) 

F - L(A ± ) / S(A)dA 
Case (ii) 

A = A + ^y 2 L"(A ) / L'(A ) +... 
i o o o 
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§4.8 DETERMINATION OF SYSTEM RESPONSE 
Knowledge of S(A) is basic to radiometry. The energy detected by a 
laboratory source of known radiance L(A) is given by equations (4.7.1) and 
(4.7.2), to which must be added a known energy owing to the detector itself. 
Thus for each pixel j: 

(4.8.1) Q - A / / / L(A) S(A) dt do> dA + Q(d) . 

T "j X 

Assume that the laboratory radiance is constant and that the pixel is uniformly 
illuminated by it. Then: 

(4.8.2) * - [Q-Q(d)]/A e T«, 

or its equivalent, is a measured and known quantity. Equation (4.8.1) becomes: 

(4.8.3) $ - J L(X)S(X)dX . 

If L(A) is perfectly monochromatic for all A » then this ideali- 
zation can be expressed as: 

(4.8.4) L( A) - 6(A-5) . 

Equation (4.8.3) is then simply: 

(4.8.5) *(5) - / S(A) (X-C)dX = S(£). 

On the other hand, with the help of (K+l) radiances: 

W’ (i ’ 

that have sufficiently different gradients, S(A) may be determined at (K+l) 
spectral points each separated by £: 

(4.8.6) A - A + ek, (k = 0,1,2, ... ,K) . 

The short and long wavelength limits are chosen so that S(A) and S(A+eK) are 


effectively zero. 
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By Simpson's Rule: 


(4.8.7) 


*1 ' <e/ 3 > X c k L i,k<WV> 

k=0 


where : 


(4.8.8) 


Let : 


1, (k - 0 ,K) 

4, (k - 1,3,. ...K-l) 

2, (k = 2,4,... ,K-2) . 


(4.8.9) \ « (e/3)c k S k . 

In matrix notation, equation (4.8.7) becomes: 


(4.8.10) [$] - [L]-m . 

Then: 

(4.8.11) m - [L]" 1 -!*] 

gives S(A) at the discrete wavelength intervals of equation (4.8.6) with the 
help of equations (4.8.8) and (4.8.9). The gradients of the radiances L(A) 
must differ sufficiently to avoid ill-conditioning. 
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V. ATMOSPHERIC EXTINCTION 
§5.1. MONO- AND HETEROCHROMATIC MAGNITUDES. 

The detected heterochromatic radiance F of equations (4.3.2) or (4.7.2) 
is expressible in terms of monochromatic radiances L(A), where A may be the 
mean or isophotal wavelength (Table 4-2). It Is desirable to express these 
quantities in terms of mono- or heterochromatic astronomical magnitudes. In 
the notation of equations (4.1.4) and (4.5.1), and with the help of equations 

(4.7.3) and (4.7.4), the detected heterochromatic radiance is: 

(5.1.1) F(S) = / L(A) S(X) dA . 

This is: 

(5.1.2a) F(B) = L(A q ) BS [1 + L”(A o )/L(\ o ) + 

in cases (i),(iii), and: 

(5.1.2b) F(B) = L(A i ) BS, 

in case (ii) [§4.2] . 

By equation (4.1.4), the detected heterochromatic illumination of solid 
angle ft is: 

(5.1.3) I(B,fi) = J F(B) du . 

In the sense of equation (4.1.7), suppose F(B,fi) is an average over solid angle 
Q, which is often simply the solid angle subtended by a pixel. Then: 

(5.1.4) I(B,fi) = fi F(B,fi) . 

By equations (5.1.2a,b), F(B,fi) is expressible in terms of radiances which are 
also averages over the pixel; thus: 

(5.1.5a) F(B ,fi) = L(A Q ,fi) SS [1 + h\s 2 L”(A o ,fi)/L(A o ,fi) + ...], 

in cases (i),(iii), and: 

(5.1.5b) F(B,fi) = L(A i ,fi) BS, 


in case ( ii) . 
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Thus heterochromatic magnitudes corresponding to the pixel solid angle 
may be calculated from measured or known quantities. By equation (4.1,8): 

(5.1.6) m(8,«) = m (0) (B) - 2.5 log 1(6, fl) . 

Since : 

(5.1.7) log(l+x) - loge (x + ...), (x«l), 
equations (5.1.4) - (5.1.6) give: 

m(g ,&) - m^°^(g) 

(5.1.8a) = -2.5 [logL(X o ,Q) + logQBS + loge V z y 2 L"(A o ,Q)/L(A o ,ft) + 

(5.1.8b) - -2.5 [logL(A i> S2) + logfigS] . 

in the various cases (i)-(iii). In this way, measured illumination is 
expressed in magnitudes. 

Evidently, it should be possible to estimate monochromatic magnitudes 
at the mean and isophotal wavelengths from the heterochromatic magnitudes of 
equations (5.1.8a, b). Let the spectral illumination incident upon the 
objective of the telescope be: 

(5.1.9) J(A,B) * J L( A ) du . 

Again in the sense of equation (4.1.7), this is: 

(5.1.10) J(A,Q) = a L(A,B) . 

Define monochromatic magnitudes at the telescope objective, and over B, to be: 

(5.1.11) m(A,Q) - m (0) (A) = -2.5 log J(X,ft), 

where n/°^(A) is the magnitude zero-point corresponding to wavelength A. Let A 

be either A or A.; then equations (5.1.10) and (5.1.11) enable L(A ,B) or 
o i o 

L(A^,Q) to be eliminated from equations (5.1.8a,b). Thus: 
m(B,Q) - m (o) (6) 

(5.1.12a) - m(A o ,^)-m (o) (A o ) - 2.5[log6S + loge ^p 2 L"(A o ,fi)/L(A o ,Q)+. . . ] , 

= m(A ,ft) - m (o) (A.) - 2.5 log6S . 


(5.1.12b) 
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§5.2. MONOCHROMATIC MAGNITUDE RELATIONS 

The relation between m(A ) and m(A.) follows directly from equations 

o 1 

(5.1.12a) and (5.1.12b): 
m(A i ) - m (0) (A i ) 

(5.2.1) = m(A o ) - m (0) U 0 ) -2.5 loge 35p 2 L"(A o )/L(A o ) + ... 

where for notational simplicity the argument of solid angle is suppressed. To 
find the corresponding relation for the effective wavelength A^, use the 
expansion AL(A) = A^LQ^) + (A-A e )[L(A e ) + A^L^A^)] + ••• to S et: 

JAL(A )S(A )dA * A e L(A e )/S(A)dA + [L(A fi ) + A g L’ (A e ) ]/ (A-A g )S(A)dA + ... 

By definitions (4.3.6) and (4.5.1), |(A-Ag)S(A)dA = (A Q -Ag)6S, so that: 
/AL(A)'S(A)dA = [X e L(A e ) + {L(A g ) + A g L' (A g )} (A 0 “A e ) IBS + ... 

Substitution of this into the definition (4.3.14a) of A , gives: 

/L(A)S(A)dA = L(-A e )[l+{L(A e )+A e L'(A e )}(A o -A e )/A e L(A e )+...]BS. 

But by equations (4.3.2) and (4.3.4), the left-hand side is just L(A.)$S. So: 

(5.2.2) L(A i ) = L(A g ) [1 + {1 + AgL' (Ag)/L(Ag)}{ (A Q /Ag)-1} + ...]. 
Conversion to magnitudes is accomplished by using equations (5.1.10) and 
(5.1.11), and the approximation (5.1.7) gives: 

m(A 1 ) - 

(5.2.3) = m( Ag) - m (0) (A e ) -2 . 51oge [ (A Q /Ag)-1 Ul+AgL’ ( Ag)/L( A g ) 1 + ... 
Nevertheless, it is desirable to express all monochromatic magnitudes 

in terms of quantities that depend only on A q . By equation (4.3.14b): 

1 — A /A = y 2 L f (A ) / A L(A ) + ««« • 
o e o o o 

Also L(A ) = L(A ) + (A -A 0 )L'(A o ) + .... and L'(Ag) = L'(A q ) + ... . Further 

reduction gives: 

m(A ) - m (o) (AJ 
e e 

= m(A Q ) - ™ (0) (A q ) 

-1 . 251oge p 2 [L"(A o )/L(A o )+2L’(A o )/A o L(A o )+2{L'(A o )/L(A o )} 2 ] . 


(5.2.4) 
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§5.3. MONOCHROMATIC EXTINCTION. 

The rule by which monochromatic radiance is diminished by the atmosphere is 

simply: 

(5.3.1) L(X,z) » L o (A) exp[-k(X)M(z)]. 

Here and henceforth, subscripts M o” denotes values outside the atmosphere; k(X) is 
the monochromatic extinction coefficient; and M(z) signifies the relative air mass 
at zenith angle z in units of the air mass at the zenith. For many practical 
purposes: 

(5.3.2) M(z) * sec z 

to an accuracy of 1% up to z * 65°; more accurate values are tabulated by 
Schoenberg (1929), and may be calculated from Bemporad’s formula: 

(5.3.3) M(z) = sec z - a^sec z -1) - a 2 (sec z -l) a - a 3 (sec z -l) 3 , 
where a^ (1*1, 2, 3) are small coefficients (Hardle 1962). Equation (5.3.3) is 
valid to z » 80° with an accuracy of perhaps 1%. 

Generalize the definition (5.1.11) to cover unextinguished monochromatic 
magnitudes : 

(5.3.4) m (A,8) - m (o) (A) = -2.5 log GL^A.fl). 

Similarly, for extinguished magnitudes: 

(5.3.5) m(A,Q,z) - m (o) (A) « -2.5 log fiL(A,fl,z) 

■ -2.5 log 8L Q (A,8)exp[-k(A)M(z) ] . 

In particular, A may be either the mean or isophotal wavelength. Subtraction and 
suppression of the argument 8 gives: 

(5.3.6) m(A,z) - m Q (A) + 2.51oge k(A)M(z) 

with 2.51oge = 1.086. This and equation (5.3.1) is the Lambert-Bouguer Law; its 
simplicity tempts photometrists to seek a similar relation for heterochro- 
matic extinction, but the wider the bandpass the more corrections are required. 
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§5.4. BANDPASS EFFECTS ON EXTINCTION. 

The heterochromatic magnitude m(3,ft) of equation (5.1.6) corresponds to 
a detected illumination over solid angle ft and bandpass 3 ; in the case of 
atmospheric extinction, it should be rewritten as: 

(5.4.1) m(3,ft,z) - m (o) (3) - -2.5 log 1(3,0, z), 

in order to signify that 1(3, ft, z) is detected when the telescope points at an 
angle z from the zenith. The definition of I must be modified to account for 
zenith angle. By equations (5.1.3) and (5.1.4), we let: 

(5.4.2) 1(3,0, z) ” J F(8,z) da - 0 F(B,z) , 

where as before we are concerned with radiances averaged over the smallest 
resolution element. Thus with the help of equation (5.1.2a) [cases (i) and 
(iii)], the extinguished illumination is: 

(5.4.3) I(g,0,z) = L(A 0 ,z) fi8S [1 +VL"U 0 ,z)/L(A o ,z) + ...], 
while in case (ii): 

(5.4.4) 1(3, ft, z) » L(A 1 »z) ftps . 

In equation (5.4.3), the approximation in square brackets involves the detected 

monochromatic radiance and its second derivative evaluated at A q . 

With the help of equation (5.3.5), the extinguished monochromatic 

radiance L(A q »z) or L(A i> z) in equations (5.4.3) or (5.4.4) may be expressed 

readily in terms of the monochromatic extinction coefficient k(A Q ) or k(A^) 

respectively. However in the case of equation (5.4.3), the heterochromatic 

illumination is also dependent on the 2 n< * derivative of the detected radiance 

L(A ,z); and this 2 nc ^ derivative will in turn entail knowledge of the 2 n< * 
o 

derivative of k(A ) as well as L M (A ). 



becomes : 
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By use of equations (5.3.1), (5.4.3), and (5.1.7), equation (5.4.1) 


(5.4.5) m(8,ft,z) - m^°^(8) 

- -2.5 log SS -2.5 log [QL(X o )] 

„ , 2 d 2 r T x -k(X )M (z)i /T > -k(X )M(z) . 

-2 . 51oge ^y 2 {h^(X^)e ° 6 ° +•** 

o 

But equation (5.3.5) can eliminate the second term in favor of extinguished 


monochromatic magnitudes at the mean wavelength. It follows that: 

(5.4.6) m(B,0,z) - m (o) (B) 


m(X ,G,z) - m (o) (X ) - 2.5 log BS 
o o 


„ „ , , d a f T ,, x -k(X )M(z) 1/t . -k(X )M( z) 

2.51oge hv 2 — 3 {L n (X n )e o' UW e ° +• 


dA ‘ l o' o 
o 


In the event that there is no air mass, M(z) = 0, and equation (5.4.6) becomes: 

(5.4.7) » o (B,8) - m (o) (B) 

- m (X ,0) - m (o) (X ) - 2.51ogBS 
o o o 

— 2.51oge ^y 2 L» o (A^)/L^(A^) + ... 

Again, the argument of solid angle may be omitted in equations (5.4.6) 

and (5.4.7), bearing in mind that it is implied for magnitudes so calculated. 

On subtracting equations (5.4.6) and (5.4.7), the magnitude zero-points and the 

bandwidth-sensitivity factor 0S fall away as they should. The heterochromatic 

extinction correction m(B,z) - ® Q (3) is then expressible in terms of the 

monochromatic extinction correction m(X ,z) - m(A ), which in turn is given by 

o o 

the equation (5.3.6). Thus: 

(5.4.8) m(B,z) - m Q (8) - 2.5 loge k(A Q )M(z) 

. 2 d 2 fT N -k(A )M(z) i t . -k(X )M(z) ^ 
-2.51oge hu 2 ^- a {L o (X n )e o' J/VV e K o' K ' +... 


o o 


+2.51oge hv 3 L"(X o )/L o CA o ) + ... 

Equation (5.4.8) reduces to (5.3.6) only for narrow-band photometry (i.e. \i 
negligibly small). 
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Inspection of equation (5.4.8) shows- that the term containing 

L"(X )/L (X ) cancels, but the term containing L'(X )/L (X ) does not. Let: 
o o o o o o o o 

(5.4.9) L ' (X )/L(X ) = lnlO dlogL (X )/dX 

o o o o oo o 

- lnlO [log L q (X+) - log L o CX o )]/[X^-X q ] + ... 
where X+>X q by the conventional use of forward differences. But by equation 
(5.4.5), the difference in the logarithms of L q at the two mean wavelengths may 
be expressed in terms of the difference in the heterochromatic magnitudes at 
these mean wavelengths; leading terms only need to be considered in this 
approximation. Thus equation (5.4.9) leads to: 

(5.4.10) 2.5 L'(X )/L(X ) 

o o o o 

- lnlO [m o (B)-m o (B + )-m (o) (6)-Hn (o) (8 + )-2.51og(B + S + /8S)]/[X^-X o ] 

+ « • • 

where : 

(5.4.11) m o (8)-m o (8 + ) = 1^(3, 8 + ) 

is a color index outside the atmosphere. The approximations (5.4.9) and 

(5.4.10) assume that the mean wavelengths X+ and X q corresponding to 
bandpasses g + and g are not too far apart. 

The required differentiations in equation (5.4.8) give the leading 
terms of the appropriate heterochromatic extinction law in terms of color 
index T , bandpass sensitivity factors gS, magnitude zero-points, and k(X Q ) and 
its derivatives k f (X Q ) and k"(A )• With the help of equation (5.4,10), the law 
is: 

(5.4.11) m(g,z) - m Q (g) 

- [k 1 (8) + k 2 (8)r o (g,8 + )]M(z) + k 3 (g)M 2 (z) + ... 
in which the new extinction coefficients are now labelled appropriately in 
broadband notation. They are: 
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(5.4.12a) 

k^B) - 2 . 51oge k(A o ) + 2.51oge ^y 2 k"(A o ) 


+ k 2 (S) [m (o) (8 + )-m (o) (8) + 2.5 log(8S/8 + S + ) ] , 

(5.4.12b) 

k 2 (B) - V a k'(A o )/(A+-A o ), 

(5.4.12c) 

- k 3 (B) - 1.25 loge [y(A o )k’(A o )] 2 . 


The leading term in the estimation of heterochromatic extinction would still 
obey Bouger f s law were it not for terms in the color index and a quadratic term 
in air mass. Conditions for neglect of these terms is derived in the next 
section. 
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§5.5. CONDITIONS FOR NEGLECT OF TERMS. 

It is instructive to estimate the various contributions to each 
extinction coefficient k (B) (i=l,2,3) in equations (5.4.12a-c). For 
illustrative purposes, choose the Johnson UBV system (Johnson 1966), whose 
system parameters are listed in the first three columns of Table 5-1. In 
addition, estimates for k(B) and its derivatives are needed. For this purpose, 
a surprisingly accurate fit to monochromatic extinction is: 

(5.5.1) k(A) = 0.24 (0.447/A) 3 . 

All wavelengths are measured in microns. Equation (5.5.1) is a quasi-Rayleigh 
scattering law, and supplies the last three columns of Table 5—1, which are 
used to derive the values in Table 5-2. 

The second and third columns of Table 5—2 account for the first two 
terms in k^B) [equation (5.4.12a)]. In all cases, the dominant term is still 
the monochromatic extinction 1.087k(A o ); the second term in equation (5.4.12a) 
is 1^ orders of magnitude less, while the last term of k^(B) is even smaller 
(since both the zero-point magnitude difference of the Johnson system, and the 
logarithm of the relative bandpass-sensitivities, are small). 

The fourth column of Table 5-2 gives the term k 2 (B) which multiplies 
the color index: 

(5.5.2) r o (B,B + ) = (U-B) q or (B-V) q 

in the forward difference convention of equation (5.4.9). The values listed 
show that the color index contribution is more than an order of magnitude 
smaller than the leading term in k^Cg), thereby warranting the label "second 
order extinction coefficient". 

The last column of Table 5-2 gives the term multiplying the square of 
the relative air mass; the smallness of k^(0) may be deceptive, since the 
product k^(g)M 2 (z) may be significant at large air masses. 
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An absolute criterion for the neglect of the color index term is that: 

(5.5.3) |k 2 (3)r Q M(z)| S e . 

For purposes of estimation, take M(z) = sec z [equation (5.3.2)]. Thus: 

(5.5.4) |rj £ e cos z / k^S), 

Table 5-3 gives the ranges of permissible color indices U-B and B-V [equation 
(5.5.2)], for e=0.01 (appropriate for 1% photometry), as a function of zenith 
angle z. It is seen that the smaller the zenith distance z, the greater the 
latitude of color indices that permit neglect of the second order color term; 
and in theory, only color indices close to zero are permitted down to the 
horizon. 

The absolute condition for neglect of the quadratic air mass correction 
is : 

(5.5.5) k 3 (B) M 2 (z) £ e, 
or: 

(5.5.6) sec z £ /e/k^S) . 

Table 5-4 gives ranges of zenith angle z over which the term may be neglected, 
for 1% or 2% photometry in U or B. Curvature in the Bouger-type linear 
relation is expected for observations that extend beyond these limits. 



Table 5-1. 

System parameters and monochromatic extinction 
terms for the Johnson UBV photometric system. 


M 2 

(pA 0 )’ 

A 

o 

k(A o ) 

k'(A Q ) 

k"(A o ) 

4.6xl0~ 4 

3.76xl0' 3 

0.350 

0.50 

-4.3 

49 

1 . 3xl0~ 3 

6.53xlO -3 

0.477 

0.24 

-1.6 

14 

1 . 6xl0~ 3 

5 . 22xlO~ 3 

0.556 

0.13 

-0.67 

4.8 


Table 5-2. 

Component contributions to heterochromatic extinction. 

A 1 . 087k( A ) 0 . 543y 2 k"(A ) k.(g) k.(B) 

O O O Z J 

U 0.350 0.54 0.012 0.020 -0.0046 

B 0.477 0.26 0.010 0.019 - 0.0018 


V 


0.556 


0.14 


0.0042 


-0.0004 






Table 5-3. 

Ranges of permissible Color Indices for 1% photometry. 


z 

0° 

10° 

O 

O 

CM 

30° 

O 

O 

<r 

O 

O 

m 

60° 

70° 

o 

O 

00 

i(U-B) o 

0.50 

0.49 

0.47 

0.43 

0.38 

0.32 

0.25 

0.17 

0.09 

±( B-V) 

0.53 

0.52 

0.49 

0.46 

0.40 

0.34 

0.26 

0.18 

0.09 

o 











Table 5-4. 

Ranges of permissible zenith distance for photometry in U and B. 
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§5.6. EFFECTIVE WAVELENGTH CORRECTION. 

The results of §5.4 et seq. pertain to case (i) (§4.2 and Table 4 — 1). 

In the event that magnitudes are required at the effective wavelength A g [case 

( iii) ] , it is seen from Table 4-1 that the previously derived magnitudes are 

unchanged, but A g must be calculated from equation (4.3.14), viz. 

A = A + p 2 L' (A )/L(A ) + . . . . 
e o o o 

This expression must now be rewritten to account for extinction. 


At the detector: 


A — A + p 2 L'(A )/L (A ) + . 

e,o o H o o o o 


(5.6.1) A e (z) = A q + p 2 L'(A o ,z)/L(A o ,z) + ... , 
whereas at the top of the atmosphere: 

(5.6.2) 

The difference between equations (5.6.1) and (5.6.2) gives the correction that 
must be applied to find the extra-atmospheric effective wavelength: 

(5.6.3) A g q — A e ( z ) = p 2 [L^(A o )/L o (A o )-L'(A o ,z)/L(A o ,z)] + ... . 

The monochromatic extinction law (5.3.1) gives the difference of the ratios: 

(5.6.4) L^(A o )/L o (A o )-L’(A o ,z)/L(A o ,z) = k'(A o ) M(z), 

and it follows that the effective wavelength correction can be written in a 
remarkably tidy form: 

(5.6.5) A 0 0 "A e (z) = p 2 k'(A Q ) M(z) . 

By equation (5.4.12b), the correction may be expressed in terms of the 
second order heterochromatic extinction coefficient: 


(5.6.6) A -A (z) = (A + -A ) k„(g) M(z) . 
e , o e o o z 

Equations (5.6.5) and (5.6.6) imply that 2-color photometry is necessary to 
derive the effective wavelength correction. The value of the correction may be 
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estimated for the Johnson UBV system with the help of Table 5-2; given that 

k («) a 0.02, and that (\ + -\ ) ^ 0.1, the correction is % 0.002 M(z), and in 
2 v p — o o 

fact need not be applied at air masses less than about 5 (z 'v 78°) in order to 
achieve an accuracy of 1%. 

The extra— atmospheric effective wavelength may be derived directly from 
equation (5.6.2), with the help of equations (5.4.10) and (5.4.11): 

(5.6.7) 

*e f o" X o 

+ 0 . 41nl0 p 2 [r 0 (8,8^)-m (0) (8)-ni (0) (B^)-2.5l0g(6V/8S)]/(X^-X 0 ) 

+ • . . « 

Xt would seem reasonable to write an analogous equation for the observed 
effective wavelength: 

(5.6.8) 

A e U> ” A 0 

+ 0.41nl0 w 2 [r(B,&^)-m (0 \B>-ni^ 0 he^)-2.51og(8 + S + /eS)]/(A^-X 0 ) 

The difference between (5.6.7) and (5.6.8) is simply: 

(5.6.9) X e 0 " X e (z) = 0 . 41nl0 p 2 [r o (6 + ,8 )~r (8 + ,6 ) ] /(X^-X Q ) , 

which at first glance seems to bear little resemblance to equation (5.6.6). 
However, with the help of the definition (5.4.11) of color index T, this 
difference becomes: 

(5.6.10) X -X (z) = 0. 41nl0 p 2 [k (g + )M(z) + ... +k ($ + )M 2 (z) 

e , o e i J 

-k 1 (g)M(z) - ... -k 3 (g)M 2 (z)]/(A^-X o ), 

in which leading terms only are retained (i.e. the second order chromatic terms 
are justifiably neglected). The terms in can be neglected too provided: 
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(5.6.11) M(z) S e)k 1 ($)/k 3 (g)y 

Equation (5.6.11) for neglect of the term in M 2 (z) differs from the criterion 
(5.5.5) inasmuch as it is a relative (not an absolute) criterion. Moreover, 
since equation (5.6.10) is already of second order, it is reasonable to take e 
to be 10%. Then from Table 5-2, M(z) must be less than about 12 in the worst 
case . 

Thus equation (5.6.10) becomes: 

(5.6.11) X -X (z) = 0 . 41nl0 p 2 M(z) [k (6 + )-k (S) ] /(X^-X ) , 

e , o e II oo 

which can be written: 

(5.6.12) X e Q -X e (z) = P 2 M(z) [k(X^)-k(X Q )]/(X^-X o ), 

with the help of equation (2.4.12a). Equations (5.6.5) and (5.6.12) are one 


and the same to first order. 
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VI. SURFACE BRIGHTNESS FROM STANDARD STARS. 

§6.1. Radiative Transfer 

If the Earth's atmosphere extinguishes light, it must also emit light. 

Let < and i be the mass extinction and emission coefficients at density p in 
v v 

the atmosphere. If the optical depth x^ of the atmosphere increases along the 
direction of an incoming light ray, then: 

(6.1.1) dx = - k p ds , 

v V 

where the pathlength s is measured outward. An Incoming ray of radiance L^ at 
frequency v is governed by the equation of radiative transfer (Goody 1964): 

dL 

(6.1.2) - j — — = L - S , 

dx v v 

v 

where : 

(6.1.3) S * j f< 

V V V 

is the source function. 

The integrating factor exp(x^) enables a formal solution to equation 

(6.1.2) to be found: 

x 

(6.1.4) L (x ) = L (O)exp(-x ) + f v S (x ! ) exp(-x +x T ) dx T . 

v v v v J v v v v v 

0 

The first term in equation (6.1.4) is simply the Bouguer-Lambert Law (5.3.1), 
according to which the incoming radiance L^(0) at the top of the atmosphere is 
attenuated by exp(-x^) by the time It reaches the telescope. When the source 
function is not negligible, the Bouguer-Lambert term must be augmented by the 
second term, which accounts for atmospheric emissivity from all points along 
the line of sight. The source function is attenuated over the optical 
thickness from the source to the telescope, and integrated over all sources to 
the top of the atmosphere. 
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If the source function were to depend solely on the radiance imposed at 
the top of the atmosphere, then equation (6.1.4) would be an integral equation 
for L^. In the present case, is completely dominated by other sources, 
whose net effect is to provide illumination to the night sky. Equation (6.1.4) 
can then be written as the sum of the attenuated radiance of a source 
observed at zenith distance z through an optical thickness t^(z), and the 
contribution from the sky "s": 

L v (z) - L v>o (*) exp[-x v (z)] + L v (z,s). 


(6.1.5) 
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§6.2. Detected Energies. 

Consider a detector that is energy cumulative, and consider each term 
of equation (6.1.5) in turn. 

(a) Celestial Sources. 

The energy arising from a source of attenuated radiance L(A,z;*) 

seen at zenith distance z is to be found from the first term of equation 
(6.1.5). Following equation (4.1.6), the attenuated starlight supplies an 
energy to solid angle oo . of pixel j through a bandpass 3 at zenith distance z 
of: 

Q.(B,oi.,z;*) = / / / / S(A) L(A,z;*) dA doj dt da . 

J 3 A T gj . S 

e 3 

Assume that the integrand is independent of areal effects; then: 

(6.2.1) Q.(B,<u.,z;*) = A / J J S(A) L(A,z;*) dA dcu dt . 

J 3 6 T “j S 

By equation (5.1.1), the detected heterochromatic radiance of the attenuated 
starlight is: 

(6.2.2) F(B,z;*) - J S(A) L(A,z;*) dA , 

6 

so that equation (6.2.1) becomes: 

(6.2.3) Q.(B,w.,z;*) = A J J F(B,z;*) doj dt . 

3 3 T a. 

3 

By equation (5.4.2), the detected illumination of the attenuated starlight is 

(6.2.4) I(B>w.»z;*) = / F(B,z;*) dcj , 

J 0) . 

J 

so that equation (6.2.3) becomes: 

(6.2.5) Q.(3,<u.,z;*) = A J I(B,w.,z;*) dt . 

J J ® t 3 


By equation (5.4.1), which is adapted from equation (4.1.8), the illumination 
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is expressed in apparent magnitudes over the solid angle as: 

-0.4[ra(8,co . ,z;*) - n/°\&)] 

( 6 . 2 . 6 ) I( 8 ,^,z;*) = 10 J 

Consequently, equation (5.2.5) is: 

-0.4[m(B,w.,z;*) “ m (0) ( B )] 

(6.2.7) Q,(B,« z;*) = A / 10 3 dt . 

J J e T 

The detected energy can be expressed in terms of magnitudes outside the 
atmosphere by use of equation (5.4.11); the argument of solid angle was 
suppressed there, and may now be reinstated, to give: 

( 6 . 2 . 8 ) Qj(Mj» z S*) 

-0.4 [m (6,w.;*)-m (o) (6)+{k 1 (B)+k (B)r(B,6 + ,cu )}M( 2 )+k.( 6 )M z (z)] 

= A J 10 ° J 12 j dT. 

6 T 

Note that the argument of m Q in equation (6.2.8) clearly implies that 

m is the unattenuated magnitude merely of that part of an image that happens 
o 

to fall in the pixel solid angle cu . Thus if the source is a star, would 
equal the unattenuated star magnitude either when the pixel solid angle is 
large enough to accomodate the attenuated (and therefore blurred) image, or if 
many pixels are used to cover the image. The former case presents no further 
complication . 

In the latter case, summation over a sufficient number of pixels, say: 

(6.2.9) j = 1,2, ...,J , 

must occur in order to ensure that all of the attenuated energy is accumulated. 
In the case of stars, their apparent sizes have been discussed by Young (1974), 
from which the value of J may be derived. Let: 

J 

(6.2.10) n = l W. , 

j-1 J 

be the solid angle, generously configured to encompass the image. As a result, 
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equations (6.2.8) and (6.2.10) lead to: 


( 6 . 2 . 11 ) 

J 

Q(6,Q,z;*) = l Q.(3,u ; Z ) 

j=l J 3 

J -0 . 4 [m (B,w.;*)-m ( ° ) (0)+{k 1 (e)+k 9 (B)r(B,6 + ,a J )} M (z)+k (6)M 2 (z)] 
_ A Y f in 0 J 1 1 J J 


A further property of equations (6.2.8) and (6.2.11) is noteworthy. 
Factors in the integrand arising from the first two terms of the exponent are 
(or may) be constant in time, but the terms in air mass M(z) are generally not 
because z is a function of time. Moreover, even at a given instant, M(z) 
varies across a pixel or a star image. Therefore make the assumptions (whose 
bounds of validity will be discussed later): 

(6.2.12a) The integration time T is sufficiently brief that a 

mean air mass can replace M(z) in equation (6.2.11). 


and : 


(6.2.12b) Air mass varies insignificantly 

across pixels and star images. 

Without further ado, equation (6.2.8) becomes: 


(6.2.13) QjCMj****) 

-0 . 4 [m (3,c,j.;*)-m (0) (B)+{k (0)+k 2 (B)r(8,8 + ,^ )}M(z)+k 3 (B)M 2 (z)] 

= A T 10 ° 3 J 

e 

while equation (6.2.11) becomes: 


(6.2.14) 


J 

Q(M,z;*) = l Q,(P. «,,*;*) 
j=l 3 3 

-0. 4 [m (B,n;*)-m (0) (B)+{k 1 (B)+k (B)r(8,8 + ,fi)}M(z)+k (B)M 2 (z)] 
= A T 10 ° 1 


e 

Equation (6.2.14) is a summation over j of the preceding equation, and 
so it must be possible to relate the magnitudes occurring in both. Taking the 
ratio in light of assumption (6.2.12b), it follows that: 
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-0.4[m (6,<u.;*)] 

(6.2.15) Q(B,8,z;*)/Q(B,« ,z;*) =10 ° J , 

in which the z-dependence is redundant given the cancellation of extinction 
corrections. In fact, a similar quotient may be formed from equation (6.2.6) 
and its analogue in ft, so equation (6.2.15) holds for apparent magnitudes as 
well. 

In reality, (6.2.15) is a flexible equation describing the relation 
between magnitudes detected by different solid angles. Only the solid angle 
arguments need retained in an obvious notation: 

-0 .4 [m(ft; *)-m( w . ; *) ] 

(6.2.16) Q(«;*)/Q j (w ;*) = 10 J . 

With the help of equations (6.2.6) and (6.2.4), the quotient is: 

-0. 4 [m( ft; *)-m(& ; *) ] 

(6.2.17) 10 J = / F(w;*)d^ / / F(w;*)do; . 

ft c o 

3 

Suppose ft comprises J pixels of size <*)y By analogy with equation (6.2.10), 
the numerator can be written: 


(6.2.18) 


/ F(w;*)d^ = I / F(tt;*)d<u . 

ft 1=1 0 ), 

j 


With a resolution of oj^ one can do no better than to let $ (*) be an average 
for the pixel. Then: 

(6.2.19) / F(<J;*)d^ = *.(*) cu . 

co 3 3 

3 

The quotient (6.2.16) becomes: 


( 6 . 2 . 20 ) 


Q( fi ;*)/Q .(“.;*) = I «. *.(*) / *>, ♦.(*) . 


j j 


j-i J J 


j j 


If a source has a uniform surface brightness, then all $j(*) T s are equal, and: 


( 6 . 2 . 21 ) 


-0 .4 [m(ft; *)-m(<u . ; *) ] 
Q(8;*)/Q.(ay*) = 10 J 

= ft/0) 

j 
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This could also have been derived from equation (5.1.6). Thus for example, if 
denotes the sky "s" of uniform surface brightness, 10 pixels are 2.5 mag 
brighter than 1 pixel. 

(b) The Sky. 

Similar arguments pertain to the derivation of energies deposited by 
the sky. By ana logy with equation (6.2.7): 

-0.4(m(B,w . ,z;s) -m^°\3)J 

(6.2.22) Q.( S , w .,z;s) = A / 10 J dt . 

J J e T 

There is no need to allow for extinction, since this is already accounted for 
by the transfer solution (6.1.5). Again invoking assumptions (6.2. 12a, b), 
equation (6.2.22) gives an equation analogous to (6.2.13): 

-0.4 [m( 6 , (*> . , z; s)-n/° ^ (S) ] 

(6.2.23) Q,(M ,z;s) = A T 10 J 

J J ^ 

If attenuated stellar energy over solid angle 8 i s recorded in addition, then 
equation (6.2.23) must be summed over all relevant pixels: 

J -0.4 [m(B ,<u . ,z ; s)-m^°^ ( 3) ] 

(6.2.24) Q(B,a,z;s) = A T I 10 3 

j=l 


(c) Detector. 

In general, the detector "d" supplies power P^(w^;d) to each pixel. 
The total energy accumulated by each pixel over integration time T is: 


(6.2.25) 


Q 4 (^;d) - / P.(^.;d) dt , 


while over solid angle ^ the integrated energy is: 


J 

(6.2.26) Q(^;d) = l J P.(*J ;d) dt . 

j=l T 3 3 


This quantity is regarded as known. 
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The energy contribution (6.2.25) from the detector (that integrates in 
bandpass 5) can be expressed in equivalent magnitudes m(oj^.;d): 

-0.4[m(w.;d)-m (o) (B)] 

(6.2.27) Q j( W j $ d) = A e T 10 J 

0 (w :d)/A T has the units of illumination as required by the definition of 

j j e 

magnitude, and is determinable under conditions of no incident light. For 
solid angle Q, the integrated dark current of the detector provides energy: 

J -0.4[m(w.;d)-m (o) (6)] 

(6.2.28) Q(Q;d) - AT J 10 2 

J-l 

which must be the same as that in equation (6.2.27) if u) ^ 35 fl there; i.e,: 

-0.4[m(Q;d)-m (o) (8)] 

(6.2.29) Q(fl;d) = A g T 10 

Equations (6.2.28) and (6.2.29) (and others) also provide the rule: 

-0.4 m(Q) J -0.4 m(oj ,) 

= I 10 J . 

j-l 


(6.2.30) 


10 
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§6.3 Formal Extinction Solutions. 


The difference: 


(6.3.1) Q.(S,w.,z;*) = Q.(3,w ,z;*+s+d) - Q (B,w ,z;s) - Q.(«.;d) , 

3 J 3 J JJ jj 

gives the attenuated energy of alone, and this in turn is a function of the 

extinction coefficients by equation (6.2.13). However, Q j( > z » s ) cannot be 
measured in the direction of The standard way around this problem is to 

measure the sky in the vicinity of and to strike a suitable average. 

Assume that: 


(6.3.2) 


pixels can be located near the source whose brightness 
equals that of the sky in the direction of the source. 


Denote such a pixel by the subscript j'; its energy is: 


(6.3.3) 


Q j ,( 6 , V ,z;s+d) = Q j ,(e, V ,z;s) + Q j ,(aj j ,;d) ■ 


By assumption (6.3.2): 
(6.3.4) ( 


but generally: 
(6.3.5) 


»z;s) = Q j( 3 ’ CU j » z; s) > 


Q j I ( I i d ) i Qj("j ;d ) • 


Equation (6.3.3) is therefore: 


(6.3.6) 


Q,t ,z;s+d) = Q (S,« ,z;s) + Q. t ( w .,;d) . 

J J J J J J 


Equations (6.3.6) and (6.3.1) give: 


(6.3.7) 


Q (B,« z;*) 


= Q.(B,oj Z ;*+s+d)-Q (B,u z;s+d)-{Q (w ;d)-Q (» ,,;d)}. 

JJ JJ JJJJ 


Quantities on the r.h.s. are now all potentially knowable; represent them by 


the difference: 


AQ jj ,(e, V w j” z; * ) 


= Q.(B,w.,z;*+s+d)-Q ,(B,6» ,z;s+d)-{Q (« ;d)-Q ,(« ,,;d)}. 

31 33 J J J J 

= (^(B.fti ,z;*). 


(6.3.8) 
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Comparison of equations (6.2.13) and (6.3.8) gives: 

+0 . 4 [m ( 8 , w . ; *)-m^°^ ( 8) ] 

(6.3.9) (l/A e T) 10 ° 3 ,Z ’* ) 

-0.4[{k 1 (3)+k„(6)r .(S,8 + ,w ) }M(z)+k (B)M 2 (z)] 

= 10 1 2 o,j j 3 

( j = 1 > 2 , . . . * J) • 

M(z), r and quantities on the l.h.s. of equation (6.3.9), are know- 
able when pixel j records energy from a source of known magnitude and color. 

The color index subscript j and argument allows temporarily for the possib- 
ility that the source is intrinsically extended; it will be dropped for sources 
that are extrinsically extended (i.e. blurred by seeing). Also, the arguments 

of Aq . make it clear that z in M(z) is appropriate to each j. Difficulty in 
JJ 

choosing pixel j' diminishes the usefulness of standard sources of large 
angular extent. Resolved standard sources are more likely to be galaxies or 
nebulae of small angular extent, and as such will be governed by the theory for 
point sources. 

For standard stars, the exponential on the r.h.s. of equation (6.3.9) 
may be taken to the r.h.s. and the equation summed over J to give: 

+0 . 4 [ m (8,^;*)-n/°^(0) ] 

(6.3.10) (1/A T) 10 ° Aq(6,Q,J1' , Z ;*) 

-0.4[{k 1 (B)+k (6)T (6,$ + ) }M(z)+k (B)M 2 (z)] 

= 10 1 z 

where summation of equation (6.3.8) has given: 

(6.3.11) AQ(B,a,fi* ,z;*) 

= Q( 6 , z; *+s+d)-Q ' (B,&* ,z;s+d)-{Q(^;d)-Q' (^' ;d) } . 

= q(B,fl,z;*) . 

m and T in equation (6.3.10) are now the unattenuated magnitude and color 
o o 

index of a star seen through relative air mass M(z). If these and Aq( Z ;*) are 
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known, the three unknown extinction coefficients on the r.h.s. require at least 
three such stars for their determination. 

Let subscripts: 


(6.3.12) i = 1 , 2 ,...,! , 


denote I standard stars at zenith distances z^. Insert subscript i in 
equations (6.3.10) and (6.3.11), and let Z^(z^) be -2.5 times the logarithm of 
the r.h.s. of equation (6.3.10): 

(6.3.13) Z 1 (z 1 ) 

+0 . 4 [ m .(M.;*)-m (o) (B)] 

=-2.5 log [(l/A e T) 10 °» 1 a Q i ( e » fl 1 ' Q, » z i ; *) J 

= - m Q 1 (B,Q i ;*) + m (0) (B) - 2.5 log [ AQ ± (B ,8* , z± ; *) /AJ] . 

Thus equation (6.3.10) becomes: 


(6.3.14) Z 1 (z i ) = {k 1 (8)+k 2 (B)r o ^(6,e + )}M(z i )+k 3 (3)M 2 (z i ), (i = 1,2,...,I). 
If M(z^) is regarded as known via equations (5.3.2) or (5.3.3), then this is a 
system of I simultaneous non-homogeneous linear algebraic equations in the 3 
unknown extinction coefficients k^($), & = 1,2,3. 

A minimum condition for the solution of equations (6.3.14) is that 
1=3; but this condition is only necessary, and not sufficient, as the follow- 
ing case demonstrates. Let 1=3 stars have the same zenith angle z. Subtract 
the I t ^ equation from those with i = 1 and 2. System (6.3.14) is then: 


(6.3.15) Z i (z) = k 2 (B) M(z) { r 0>i (M+) - r 0>3 (M + >} , i = 1,2. 

Only the secondary extinction coefficient can be found, and then only if the 
color indices of stars i * 1,2 are sufficiently different from the third; the 
extinction coefficient k^ which accounts for most of the extinction, goes un- 
determined. This leads to the important result: a necessary and sufficient 
minimum condition for the solution of the 3 extinction coefficients is that 
there be 3 standard stars with different zenith angles and color indices. 
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In practice, it is necessary to overdetermine the system (6.3.14) by 
use of I >> 3 standard stars having a large range of colors, and to proceed by 
use of least squares. Since MCz^) is common to all terms on the r.h.s. of 
equation (6.3.14), and is assumed to be known, the system to be solved may be 
written as: 

(6.3.16) Z 1 (z i )/M(z 1 ) = k 1 (e)+k 2 (6)r o>1 (6,g + )+k 3 (B)M(z i ), (i = 1,2,...,!). 
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§6.4. Zenith Angle and Color Index Ranges. 

A result of §6.3 is that a range of z is necessary for the solution of 
of the primary extinction coefficient k^ in system (6.3.14), whereas a range of 
r is necessary for the solution of the secondary color coefficient k 2 • More- 
over, the tertiary coefficient waits to assert itself at large z. Such effects 
represent competing claims to significance, and thus to numerical accuracy. To 
guage the situation in a preliminary way, consider the following arguments. 

The range of color index affects only k 2> whereas the range in air mass 
M(z) is common to all terms in equation (6.3.14). Thus, let: 

(6.4.1) z ± - z ± «z ± , 

where z is a zenith angle somewhere in the middle of the range of the I stars, 
and 2 1 <5 Z | < F.O.V. (the field of view of the detector). To evaluate the sizes 
of the competing claims, it suffices to take M(z) = sec z. Then: 

(6.4.2) M(z i ) = sec(z±5 Zi ), 

which can be written: 

(6.4.3) M(z i ) = M(z) {l ± f i (z,6z i )} . 

To a first approximation, let: 

(6.4.4) M(z) = sec z , 

so that the function f is exactly: 

(6.4.5) f^(z,&z ± ) = (cos - tan z sin 1 - 1. 

Values of f are listed in Table 6-1. In particular: 

(6.4.6) f.(z,0) = 0, 

and : 

(6.4.7) f (0,6z ) = sec 6z - 1 . 

Choose z to be the zenith distance of the 1^ star, so that = f (z,0) = 0. 

Applying equation (6.3.14) gives: 



Table 6-1 
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(6.4.3) Z.(z.) - Z I (z I ) 

= ±k 1 M(z)f 1 + k 2 M(z)[(i±f.)r o>i - r o>I ] + k 3 M 2 (z)(±2f i +fp. 

When = 0, the case pertaining to equation (6.3.15) is recovered 
(i.e. the primary coefficient can be found only if there is a range in z). But 
| f i | must then not be too large lest it be swamped by the tertiary contribu- 
tions. Generally it would be desirable to have: 

(6.4.9) k x 2|k 3 M(z)] | 2±f 1 |. 

By Table 5—2 , J ^ is two orders of magnitude less than in the Johnson system, 
so | 2±f^ | £ 100 cosz or: 

(6.4.10) 2 - 100 cos z < |f^| < 2 + 100 cos z . 

The more stringent inequality on the l.h.s (for stars closer than average to 
the horizon) gives the absolute limit of no permissible range in z: cosz = 0.02 

for an elevation angle of l°10 f . Only narrow ranges of z are permissible if 

bandpass effects beyond y 2 are neglected, whereas near the zenith, large 
ranges are allowed [and are desirable in fact, according to equation (6.3.15)]. 

Competing requirements on ranges of z and color are evident in the sec- 
ond term of equation (6.4.8). Consider the conditions on the first and second 
terms of (6.4.8), viz. that (i) the second term doesn’t dominate the first, and 
conversely (ii) the first doesn’t dominate the second. For want of a number, 
write these respectively as: 

(6.4.H) kjf.l 2 |k 2 [(i±f i )r 1 -r I ]| fckjfj/ioo. 

By Table 5-2, k^ is an order of magnitude less than k^ , so that: 

(6.4.12) 10|f.| £ U r i - r I ) ±f i r i l s l f ± 1/10 . 

Consider two cases: 

(a) if |f^| is small compared to | | then: 

iolf.1 2 I r 1 -r i | S IfJ/io . 


(6.4.13) 
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The case of equation (6.3.15) is implicit in the first inequality, which 
states also that a small range of z restricts the simultaneous solution for the 
primary and color coefficients. The second inequality states that the range of 
color index differences must exceed a certain fraction of | f ^ j in order for the 
second term to be determinable to a given accuracy. 

(b) if |f | is large compared to Ir^-T^l, then: 

(6.4.14) 10 2 |r | > 1/10 , 

which is reasonable. 
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§6.5. Apparent Magnitudes. 

To describe the appearance of a pixel in magnitudes m( g , z ; *+s) , it 
is necessary to subtract detector energies QU^d) from total energies 
Q( 8 > w • > z > *+s+d ) 5 i • e . ■ 

J / \ 

-0.4[m(8,<u .,z;*+s)-t/ 0 ; (g)] 

(6.5.1) Q(g, Wj , 2 ;*+s+d)-Q( Wj ;d) = A e T 10 3 

With the help of equation (6.3.8), this can be expressed in terms of AQ(*)= 

-0.4[m(g, w . ,z;*+s)-n/°\g)] 

(6.5.2) A g T 10 J 

= AQ(e,w.,u,. ,z;*+s+d) + Q'( 8 , 0 ),, ,z;s+d) - Q'^.id). 

J J J J 

The appearance of the source follows from summation over J. 
-0.4[m(£,Q,z;*+s)-m^ O \$) ] 

(6.5.3) A g T 10 

= Q(B>ft> z ;* +s+ d) ■ Q(fl;d) 

= AQ(8t0>n' ,z;*+s+d) + QHe.fl' ,z;s+d) - Q'(fl';d). 

In particular , if ,, * M contributes nothing to the differences In 
equations ( 6 . 5 . 1 ) and ( 6 . 5 . 2 ), then the apparent magnitude of the sky follows 
from: 

-0 . 4 [m( 3 ,oj • > z , s ) m ( 3 ) ] 

(6.5.4) A e T 10 3 

= Q ' ( 8 , w . , , z ; s+d ) - Q'(aj.,;d). 

J 

More reliable statistics for sky brightness in the vicinity of a star result 
from summation over J: 

-0.4[m(g,Q,z;s)-m^°^(B)] 

(6.5.5) A g T 10 

= Q' ( 8 , £2' ,z;s+d) - Q'(fl';d). 

Equation (6.2.30) may be used to reduce sky brightness to convenient units 
(e.g. mag /arc sec 2 ) . 
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§6*6 Derived Magnitudes. 

In equation (6.3.13), let subscript "u" connote a source whose un- 
attenuated magnitude is unknown. Then by writing the color index explicitly 
according to equation (5.5.2), the unattenuated magnitude in bandpass g is: 

(6.6.1) ra o u (8,n,*u) = n/°^($) - 2.5 log [ AQ ( & , T ,Z u ;*u)/A e T] 

- k 1 (g)M(z u ) - k 2 (8)M(z u )[ni 0>u (6)-ni 0>u (6 + )] 

- k 3 (g)M 2 (z u ) . 

The extinction coefficients can be determined from I>>3 standard stars, and 
only the forward component in bandpass g of the color index [m Q u ^B)” m 0 U (B )] 
remains to be determined. To get ^(g ) requires a second solution in g 
from a second exposure of time T and zenith angles z ; thus: 

(6.6.2) m (B + ,a,*u) = m (o) (B + ) " 2.5 log [AQ(3 + ,Q,«' ,z^;*u)/A T + ] 

- k 1 ( g + )M( Z ^) - k 2 (e + )M(z^)[m 0 jU (B + )-m o>u (B^)] 

- k 3 (B + )M 2 (z^) . 

The need for m Q U (B + ) in equation (6.6.1) arises only in the second 
order color coefficient, and thus only first order terms in equation (6.6.2) 
are necessary. Thus replace m^ ^(g + ) i n equation (6.6.1) by: 

(6.6.3) m Q u (B + ) = m (o) (g + ) - 2.5 log [AQ(B + ,G,Q' , z*; *u) /A g T + ] 

- k 3 (g + )M 2 (z^) » 

where the tertiary coefficient term is retained lest it be significant. Thus 
m o u (g,ft,*u) can be determined. 

To determine the unattenuated color index of the source, the difference 
between equations (6.6.1) and (6.6.2) is necessary; but to get m Q ^(g ,ft,*u), 

I | | 

it is necessary to know the next forward magnitude m Q U (B )• This may be 
found in the same way from a third solution in bandpass g . Thus in the 
Johnson U,B,V system for example, unattenuated magnitudes (U,B) and color 
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indices (U-B,B-V) can be found. Then magnitude: 

(6.6.4) V = B - (B-V) 

follows . 

Instead of using the forward convention of equation (5.4.9), backward 
differences could also be used, in which case the analysis of equations (5.4.9) 
to (5.4.12) must be reworked. Then the backward difference equivalent of 
equation (6.6.4) would be: U = B + (U-B). It is noteworthy that the hetero- 
chromatic coefficients k^(8) and k^(B) would then be slightly different, as 
seen from equations (5.4. 12a, b). Presumably, the magnitudes obtained by for- 
ward differences would be the same as those found by backward differences to 
the accuracy retained in the analysis. 

Equations (6.6.1) to (6.6.3) obtain equally for pixel solid angles oj y 
Thus unattenuated surface brightnesses and colors of extended sources can be 
obtained for each pixel, and reduced to convenient units (mags/arcsec 2 ) with 
the help of equation (6.2.21). 
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VII. DETECTOR CALIBRATION . 

§7.1 Analog Devices. 

Energies Q appearing from equations (4.1.6) to (6.6.3) are required in 
digital form in order for solutions to be found, yet often are recorded by 
analog devices. Such devices are usually not precisely linear; an extreme 
example is the photographic emulsion, though in both CCD and photographic 
photometry the approach to saturation manifests itself as a curvature in the 
response . 

In general, the energy deposited on pixel j may be expressed as a 
function of a dimensionless parameter D that measures the analog output: 

(7.1.1) Q d - Q j( D.) . 

In the case of photography, D would be the specular photographic density. 
Following the work of Zou, Chen, and Peterson (1981), let the function Qj(D^.) 
be modeled by a polynomial series: 


N 


(7.1.2) 




where parameter N is assigned and p is an optimization parameter that minimizes 

the discrepancy between observed and derived magnitudes of the standard stars 

( o ) , 

for the assigned N. Note specifically that the coefficients C^ in equation 

(7.1.2) have units of energy. In the case of photography, the polynomial P may 
be the standard rectification function of Baker (1925) and Sampson (1925) 
raised to the power (np): 


D . np 

(7.1.3) P n . (D ; p ) = [10 J - 1] . 

Wherever energy Q has appeared in the foregoing sections, it may be 
replaced by the development (7.1.2). The relevant notation signifies the 
conditions that obtain; thus for example, subscript j and argument "s+d" imply 
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that Dj(s+d) is the analog output for pixel j for sky plus detector* To re- 
lieve notational overload, the bandpass arguments will be dropped * 

Thus from equation (6.3.8) it follows that: 

(7.1.4) A Qj j* ,z; *) 


N 


* ^ C n 0) r p n>j {V* +s+d >;p} - p n>j .{D r (s+d);p} 

- P n>j {D.(d);p} + P n> .,{D jt (d);p}]. 
However by equation (6.3.9), this is: 


(o) 


(7.1.5) A T 10 
e 

N 


-0.4 [m (oj . ;*) - itT ' + k_M(z) + k 0 M(z)T + k„M 2 (z) 

O j i ZOO 


■ l C n 0> ,P n,jl D j < * +s+d);p l ' p „,j-l D j .< s+d >;Pl 

- + p „ >J - { V d)ipl1 

Since coefficients are to be determined, they might as well be 


C = (1/A T) 10 _0 * 4m C (o) 
n e n 


replaced by: 

(7,1.6) 

Thus the effective collecting area is absorbed into determinable coefficients; 
its value need not be known, but what is required is that: 


A is the same wherever a 

0 

solution is to be applied. 


(7.1.7) 

Note that this limitation imposes a further constraint on the range of zenith 
angles (cf. §6.4). 

Equation (7.1.5) becomes therefore: 


(7.1.8) 10 


-0.4 [m o (a> ;*) + k^z) + k 2 M(z)T o + k 3 M 2 (z)] 


N 


= I C n f P n , j { D j (* +s+d );p} - P n> j i { D j * ( s+d) » p} 

- P n>j {D j (d);p} + P n)jI {D jt (d);p}]. 
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Summing over j yields an expression for the known magnitude m . (*) of star i: 

0,1 

-0.4 [m 0j .(*) + k 1 M(z 1 ) + V 1(Z i )r o,i + k 3 M2(Z ! )] 


(7.1.9) 


10 


■J^n j j„ 1 [P n,:,i lD 3,i < * +S+d,;,>l - P n,j',i ,D J'.i (S+d):pl 


(i = 1,2,...,!). 


This is a nonlinear system of equations that permits solution for the three 
extinction ceofficients and the N calibration coefficients. 

Equation (7.1.9) is the denoument of the foregoing analysis. Amongst 
other things, it shows that ground-based detector calibration cannot be accom- 
plished by use of standard stars unless extinction coefficients are simultan- 
eously part of the solution. This principle has previously been asserted 
(Usher 1986), but is here spelled out in complete detail, along with implicit 
assumptions. In addition, the apparent disappearance of the magnitude zero- 
point in the work of Zou et al., is clarified, as is the role of the effective 
collecting area. 

Moreover, the sky brightness will now emerge as a derived quantity, 
whose value also cannot be found without a concomitant solution for extinction; 
again, this principle is in accord with the equation of radiative transfer and 
its solution (6.1.5). From equation (6.5.4): 

-0.4m(^ . . ,z; s) N 
(7.1.10) 10 J = l C 

n=l 


[P {d ,(s+d);p} - P n {d (d) ;p} ] , 

11 11 J J J 11 J J J 


or from equation (6.5.5): 


-0.4m(ft' ,z;s) N J 

10 - I C I [P n . , {d j ' ( s+d) ; p } - P {d (d);p}]. 

n-1 j'-l ,J J 


(7.1.11) 
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The appearance of pixels (i.e. with detector noise subtracted) is found 


from equation (6.5.1): 

-0 . 4m(oj . , z ; s ) N 

(7.1.12) 10 J A = I C [P {D (*+s+d);p} - P | 

q U n >JJ 11 >J 

n= 1 

The apparent (attenuated) magnitude of a star i (i.e. with sky 
noise subtracted) follows from equations (7.1.9) and (5.4.11): 


(7.1.13) 10 


-0.4 m. (z ;*) 
i i 


= I C l [P {v ( *+s+d ) ; p } - P , ± {D 

n-1 n j,j'=l n * 3,1 J> ,J ’ J 

- P . .{D. . (d ) ; p } + P ,{D., 
n, J , i 3 , i v ' ,y> n, j' , i 1 j' 


D Cd>;p}l . 
and detector 

(s+d);p} 

y 1 

,(d);p}]. 

y ^ 
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§7.2 Perturbation Iteration. 

Given the complexity of non-linear multivariate analysis, it seems 
reasonable to attempt a solution to equation (7.1.9) by linearization and 
iteration, using Newton's method. Let: 


(7.2.1) 

aP l,n (p) ", l P „,J,i( D j> l ( * +S+,i>;P t ' P .,J',il D j',i (SW);P l 

J > J 1 

- P n> J,l( D J,i <d,;Pl + P n,j',i( D JM (d);P l 1 ' 

so that equation (7.1.9) is (with self-evident notational economy): 

-0.4fm.+k, M.+k.M. r .+k.M? ] N 

„ 1 1 lx 2 1 i 3i J 


(7.2.2) 10 


I AP i n 
n=l 1,n 


(p) C R , (i = 1,2, ...I). 


Perturbation of both sets of coefficients gives: 

(7.2.3a) k £ = k< c) + 6k £ , (Z = 1,2,3), 

(7.2.3b) C = C (c) + 6C , (n = 1,2, ...,N), 

n n n 

where superscript (c) signifies a value initially guessed or subsequently 

( c ) 

corrected during the iteration process. Values for k can be estimated from 

C c ) 

Table 5-2 with fair accuracy; corresponding values of can be found by 

solving equation (7.2.2), which is for this purpose a linear regression 

problem. Corrections 6k rt and 5C are therefore small, so that to first order: 

Z n 

(7.2.4) 10 X = 1 + ( lnlO) x + ... . 

Corrections 6k. and 6C will be governed by the following linearized 
Z n 

system, derived by substituting equations (7.2.3a,b) into equation (7,2.2), 
using equation (7.2.4), and retaining first order terms: 


(7.2.5) 


0 . 41nl0 M , ( 6k. +T , 5k +M. 6k^ )+10 
i 1 i 2 l 3 


0.4 [m +k^ C ^M +k^ C ^r.M +k^ C ^M? ] N 

11 I ^ 11J 1 C ATI /Nro 

• 2 4P l,n <p) 6C „ 
n= 1 


0.4[m i +k^ c) M i +k^ c) r i M 1 +k^ c) Mp N 


' l AP i 
n=l i ’ n 


(P) c 


(c) 


1 - 10 
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Equation (7.2.5) is now an overdetermined linear system of the form: 

(7.2.6) A X = Y , 

whose solution is (Bates and Watts 1988): 

(7.2.7) X = (A T A) _1 A T Y . 

"T" signifies a transposed matrix. Specifically: 

(7.2.8) X = 


(7.2.9) 


A = 


6ki 6k^ 

5C, 6C 

• • « « 

6C | T , 



1 2 


N 1 


0. 41nl0 M 1 .. 

0.41nl0 


AP u (p) ... 

ap in 

0 . 41nl0 M 2 .. 

0.41nl0 

M 2 

ap 21 (p) ... 

AP 2N 

* • < • 

0. 41nl0 Mj .. 

• • 

0 . 41nl0 


• « • 

AP n (p) ... 

ap in 


and : 


(7.2.10) 


Y = 


1-10 


0.4[ii 1 +k| c y+k^rj +k^ c) Ml ] N , v 

11 12 113 y ap, ( P )c (c) 

L In n 


n=l 


0.4[m +k (c) M +k^ c) r M +k^ c) Mn N . . 

2 1 2 2 2 2 J • y AP„ (p)(r ' 

L , 2n r n 
n=l 


1-10 


0.4[m T +kJ c) M T +4 c) r M +k* c) M?] N , . 

1-10 11 12 113 '• l 4P In (P)c' C> 


n=l 


Iterative solutions progress to the point where vector X gives no 
improvement in the corrected coefficients, whereupon successive completed 
iterations give the residual sum-of-squares function of p: 


(7.2.11) 


S(p) = l Y. 2 (p) . 


i=l 


Comparison of equations (7.2.10) and (7.2.2) shows that a best solution to the 
nonlinear regression is achieved when S(p) is a minimum. Thus p is determined 
for a given N and used in subsequent calculations. 
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§7,3 Air Mass as Unknown, 

In §7.1 and §7.2, it is assumed that the relative air mass M(z) is a 
known function (5.3.2) or (5.3.3) of secant z. However the secant law is a 
consequence of the assumption of plane parallelism, and may start to fail at 
larger zenith angles (§5.3), even though the law is still a fairly good 
approximation . 

The actual dependence M(z) may be found as part of the solution of the 
nonlinear least squares problem (7.1.9), by perturbing M(z) as well as the ex- 
tinction and polynomial coefficients; thus in addition to equations (7.2.3a,b), 
let : 

(7.3.1) M i (z 1 ) = M^ c) (z i ) + 5M jL (z i ), (i-1,2 1), 

where the initial guess might be found from equation (5.3.2). Regarding 

Mi( z i) as an unknown for every star i implies that there will always be (N+3) 
more unknowns than there are standard stars for their determination. In 
practice, corrections will be virtually the same within observational error 
for stars within a particular range of zenith angle; this range might even en- 
compass the entire field of view. Let there be: 

(7.3.2) m = 1,2, ...,M 

ranges of zenith angle for which a solution is sought, and: 


(7.3.3) i = 1,2,... ,I m 

stars In each range, so that the sum total of available standards is still: 


(7.3.4) ! = 1 + I + . . . + ! m . 

The statistical analog of Newton* s Method can then be applied to each band of 
zenith angle. 

Retention of first order perturbed terms gives the generalization of 


equation (7.2.5): 
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(7.3.5) 


0.41nlO[ 6 k 1 «< c) +6 k 2 r i «< C>+ « k 3 { f ’m‘ :) | ,+4M «!^ C)+r i k 2 C>+2k 3‘ :) ) 1 


+ 10 


0.4[m.+kJ c) M (c) +k$ c) r M (c) +k^ c) {M (c) l 2 ] N 
L il m 2 I m JimJ 

• l 4P i,n (p) SC n 
n=l 


0.4[ra,+k5 c) H (c) +k^ c) r,M (c) +k^ c) {M (c) } 2 ] N 
L i 1 m 2 i m 3 i m J ' 


1-10 


• l AP l,n < >’ ) C n 

n=l 


(i = 1,2,. ..,1; m = 1,2 M). 
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§7.4 Role of Spot Sensitometry . 

Spot sensitometry is sometimes used in the more intractable case of 
photography. Portions of the field are exposed to sources of light "S" of: 

(7.4.1) m = 1,2, 

different brightnesses. To forestall reciprocity failure, the duration of the 
sensitome tr ic exposures must equal the integration time T; and to avoid contam- 
ination, the calibration spots must be placed on a shielded portion of the 
field. 

In the case of the UKST (United Kingdom Schmidt Telescope; Tritton 
1983), the illumination is automatically switched on and off as the telescope 
shutter opens and closes; thus calibration spots are inserted under environ- 
mental conditions identical to those for the main field. Illumination is pro- 
vided both (a) by projectors powered by tungsten lamps that provide seven steps 
arranged sequentially on two edges of the field, and (b) by a projector of the 
Kitt Peak design powered by a quartz-halogen lamp with a color correction fil- 
ter that provides sixteen steps arranged in a square. The intensity of illum- 
ination may be varied to suit the exposure at hand. 

t h 

Following §6.3, let the energy deposited in the m exposure be the 
difference between the energy recorded and that contributed by the detector, 
viz . : 

(7.4.2) QfS) = Q(S +d) - 0(d) . 

mm mm m 

The notation of equation (7.4.2) is simplified by omission of the bandpass 
argument, the solid angle of the exposure (inasmuch as averages over large 
numbers of pixels j = 1,2,...,J can be made), and the primes on j (since 
detector contributions can only be assessed from pixels suitably near the 
spots). In the case of photography, energies on the r.h.s. of equation (7.4.2) 
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may be rewritten with the help of equation (7.1.2); i.e.: 


(7.4.3) 


N 


q (s ) - y 

mm L n 
n=l 


c(0 ’ 


P „m( D(d; P>H ■ 


where the extra polynomial subscript m is inserted without confusion to signify 

the spot under consideration (i.e. P still signifies P , but for spot m) . 

nm n 

Equation (7.4.3) can be put to a variety of uses. Energies may be 
known (a) absolutely, or (b) relatively; and in either case, that information 
may or may not be used in conjunction with calibration from standard stars. 

The more common (indeed perhaps the only) real case is that of relative cal- 
ibration, but the analysis proceeds naturally from the former. 

(a) Absolute Calibration. 

If Q is known absolutely, then equation comprises M equations in N 
m 


unknown coefficients, which by equation (7.1.6) are: 

.(O) 

(7.4.4) 


^(o) ^ ™ lrt 0.4m 

C v ' = C A T 10 
n n e 


Substitution into equation (7.4.3) reveals the connection between the absolute 
energy of calibration and the magnitude zero-point [e.g. equation (6.3.9)]; 
thus it is reduced to dimensionless form: 


<7 - 4 - 5) VV 10 '°' 4m <1/A e T> ■ l C n [ p „ n ( D ( S ” td ;P>) - p „„{ D < d -I > >}J • 

n=l 

If Q is known, then equation (7.4.5) provides M equations In N unknowns C . 
m n 

Since N might be perhaps 3 or 4 at most, the linear system would be overdeter- 
mined for the 7 to 16 spots available on UKST plates. 

(b) Relative Calibration. 

More reliable information is available from the ratios of supplied 
energies. Take ratios conveniently (with respect to Q^, say) and let: 


R m< s > S m> 
mM m M 


W'W 


(7.4.6) 
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Then: 


N 


y c [P { D ( S +d ; p ) } - P {D(d;p)}] 
L n nm 1 m 1 nm 1 } 


n=l 


(7.4.7a) R M (S ,S M ) = 
v y mM m M 


I C „ [ P n„! D ( S M +d; P ) l - P „Ml D < d -P>)l 


’ ( m= 1 , 2 , . . . , M~ 1 ) 


n=l 


nM 1 


or on rearranging terms (and suppressing argument D for simplicity): 


(7.4.7b) 


1 C [R M {P M ( S +d)-P (d) } - { P ( S +d ) -P ( d ) 

L n L mM 1 nM M nM J 1 nm m nm 


n=l 


= 0 . 


Since relative calibration is sought, only ratios of coefficients need to be 
found. Let: 

(7.4.8) 

Equations (7. 4.7a, b) are then: 

N-l 


K „ = C /C = C (0) /C< 0) , (n=l , 2 , . . . ,N-1) . 
nN n N n N 


y K „ [R Jp M (S M +d)-P M (d)} - {p n (S +d)-P (d)}] 

L nN 1 mM 1 nM M nM i 1 nm m nm 


n=l 


■ ( P N„'V d) - P Nm< d >l - R rt .! P ™ <S M +d) - P NM (d) l - 

(m * 1 , 2 , . . . , M- 1 ) . 

As noted above, relative spot sensitometry may be used in conjunction 
with calibration from standard stars, or independently. The two cases are re- 
presented by equations (7.4.7b) and (7.4.9) respectively. 

When both standard stars and spot sensitometry are available, equation 
(7.4.7b) comprises an additional M-l equations for simultaneous use with system 
(7.1.9); or it may be inserted simultaneously into the perturbed system (7.2.6) 
as part of the iteration procedure. Alternatively, equation (7.4.9) is an 
overdetermined linear system in (N-l) unknown relative calibration coefficients 
whose solution gives the shape but not the zero-point of the calibration curve. 
Relative surface brightnesses may be obtained, but rough absolute brightnesses 
may be inferred even if there is only one standard star in the field. 
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§7.5 Calibration from Space. 

In the absence of atmospheric emission and extinction, equation (7.1.9) 
still holds, except that terms must be reinterpreted. 

m (*) is the magnitude of source i attenuated only by the inter- 

o , i 

stellar medium; [corrections for interstellar extinction are not considered 
here. ] 

M(z) is the relative interplanetary mass at angle z measured toward the 
ecliptic. Extinction coefficients k^ k 2> and k 3> pertain now to the interplanetary 
medium. The integrated source function of equation (6.1.5) manifests itself as 
contributions "s" to the analog outputs D. 
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§7.6 Reciprocity Failure 

Reciprocity failure (RF) refers to the failure of a detector to main- 
tain a reciprocal relation between illumination I and integration time T in the 
definition of Exposure 1 : 

(7.6.1) E = IT . 

A given exposure obtains as long as I and T are inversely proportional to one 
another, but the detector might not respond identically in all cases. Photo- 
graphic emulsions are particularly susceptible. 

Wide-angle cameras employ photographic emulsions as the detector of 
choice in the study of large scale phenomena. In the case of comets with high 
proper motion, the existence of morphologically interesting structures can best 
be seen when the camera tracks the comet. Unfortunately, if photometric data 
are to be extracted as well, standard stars will appear trailed relative to the 
comet image. In this case, the detector is exposed to various levels of light 
for different periods of time, so RF may manifest itself in different responses 
to high I for short T than to reciprocally low I for long T. This section com- 
prises (i) the description of the problem of RF stripped of extraneous complex- 
ity, followed by ( i i ) a scheme to overcome the difficulties. 

On suppressing the argument of bandwidth 3, equations (6.3.8), (7.1.4) 
and (7.1.5) give: 

Qj( aJ j ,z; *) 

-0 . 4 [m («.;*) - m (o) + k^MCz) + k r M(z) + k.M 2 (z)] 

. m ~ o i 1 z o o 

= A T 10 J 

e 

cj; 0) [P n ,{d (*+s+d);p} - P • , f D , ( s+d ) ; p } 
n n , j j “jJJ 

- P {d (d) ; p} + P {d (d);p}] . 

n > J J n 5 J J 


(7.6.2) 
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To reveal the essence of RF, assume no detector noise, no atmospheric extinc- 
tion or emission, and a single pixel solid angle. Moreover, it suffices that 
the detector response be represented by a function P[D(*)]. For a star i , 
equation (7.6.2) is: 

(7.6.3) Q(* i ) = P[D(*.)] = A g T 10 - 0 - 4 [ m < . 


In terms of illumination: 


—0 4 rmf * ^ — rn 

(7.6.4) I(* ± ) = Q(* i )/A e T = P[D(* i )]/A e T =10 1 ± J 


(o). 


In terms of exposure: 

(7.6.5) E(* a ) = I(* i )T = Q(* i )/A e = P[D(* i )]/A e = T Hf 0 * 4 ^ V"® 1 

Perform a thought experiment in which two point sources are detected, 
one trailed and the other not. Assume that the trailed source results from 
faithful guiding on the untrailed source, but that it has zero albedo except 
for a spot that is a specular reflector. The trailed source "t (a rotating 
black asteroid with a mirror?) is detected as a series of starlike images, 
whereas the guided source "g" (a fixed star tracked at the sidereal rate?) is 
detected as a single image. The exposure time T of a single image of the 
trailed source is therefore shorter than the exposure time T of the guided 

O 

source: 


(7.6.6) T t ^ T g . 

Suppose next that sources "t" and "g" have the same exposure E, so that 
source ”t" would have a brighter image if it were not trailed; i.e.: 

(7.6.7) m(* t ) < m(* ) . 

The same energy is detected in the brighter source in the shorter time, as is 
detected in the fainter source in the longer time. Then if there were perfect 
reciprocity between I and T: 

— 0 Am ( A \ 

E(* t ) K* t ) T t Q( V T t 10 * t P t [D(* t )] 

«v " I( V T s ’ «V " t 8 


(7.6.8) 


1 . 
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The first four ratios are all unity by agreement and by definition of physical 
quantities, but the last ratio is unity only if the subscripts make no differ- 
ence; i.e. in the absence of RF : 

(7.6.9) P t [D(* t )] = P g [D(* g ) ] 

In other words, the quantity D for the brighter image "t M detected over a 
shorter time would equal D for the fainter image "g" detected over a longer 
time, and the perfect detector would have a common response P(D) under all 
circumstances . 

On the other hand, if RF exists, then: 

(7.6.10) P t [D(* t )] * P g [D(* g ) ] . 

In the experiment at hand, the problem can be handled only if a way is found to 
relate P to . This translates into different polynomial summations for the 
r.h.s. of equation (7.6.2), except that with atmospheric absorption and 
emission, the sky brightness and detector noise components are not subject to 
RF to a first approximation. 

The general problem is even more difficult owing to the fact that a 
trailed star image, blurred by seeing, is subject also to the Fog and Inter- 
mittency Effects. In the former effect, the leading parts of the trailed image 
pre-expose and pre-sensitize the detector, while in the latter effect, parts of 
the trail are intermittently exposed to equal brightnesses of the trailing 
image (except in the extreme lateral edges). Couple these and other likely 
effects, with the fact that RF is itself highly variable, and it is clear that 
photometric programs on faint structures in high proper motion objects must be 
designed with paramount emphasis on photometric control. 

It Is not sufficient merely to apply equal but oppositely directed 
trailing to both standard stars and comet by tracking at a rate equal to the 
sidereal rate plus half the proper motion rate. Not only will a great deal of 
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morphological information be lost or difficult to retrieve; but the photometry 
of a trailed point source will not be the same as the photometry of a diffuse 
source, because the point source will trail onto parts of the detector pre- 
exposed merely to the sky, whereas parts of the diffuse source will trail onto 
other parts of itself, as well as the sky. 

It Is equally unsatisfactory to rely solely on sensitometry spots, for 
while these are certainly helpful in acquiring relative brightnesses of the 
guided source, they do not necessarily provide the zero-point of the photometry 
in a reliable way. Moreover, if they are available at all, sensitometry spots 
are often applied under environmental conditions different than those of the 
main exposure. 

We suggest that there is only one genre of solution to the problem, 
with two chief variants. The underlying principle is well-known, viz. that the 
existence of strict controls and the use of standard stars Is the only way to 
do absolute photometry. The variants divide along the following lines, accord- 
ing as sensitometry spots are (a) utilized or (b) unavailable. 

(a) 

If sensitometry spots are available, then the goal is to use standard stars 
to calibrate them! Standard stars are therefore the indirect absolute cali- 
brators, and the best ones at that. A series of exposures of standard star 
fields would be made, and spots inserted, under variable environments but with 
identical integration times. Ideally, the laboratory conditions must be equal 
to the outside environment, and spot intensities must be strictly reproducable . 
All emulsions would be from the same batch, be developed together, and other- 
wise be subject to egalitarian treatment throughout. Data reduction would 
follow the lines of §7.4. There would be no RF problem, and indeed no need for 
standard stars to be in the field of the comet. [If there are standards in the 
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field, so much the better, for then they can be used to develop a library of RF 
solutions (7.6.10) as a function of exposure time, incident radiance, and 
detector type.] Several series through the night will monitor possible changes 
in extinction and sky brightness. If the opportunity window is limited, then 
the comet photometry should at least be bracketed by spot calibration. RF can 
be beaten by use of a detector whose unsusceptibility can be demonstrated via 
(a) above. Failing this, there is (we suggest) only one recourse. 

(b) 

If sensitometry spots are unavailable, a scheme must be devised to avoid the 
trailing of the calibrating stars. Just as sensitometry spots are inserted on- 
to the main exposure in a separate procedure, so standard stars (which now take 
the place of the spots) must be inserted separately. The basic maneuver is to 
allocate a potion of the detector for guided standard stars at the same mean 
zenith distance as the guided comet, which is exposed separately onto the rest 
of the detector. This can be accomplished by masking the relevant part of the 
detector during each exposure. Each exposure has identical integration times, 
and is made as soon as possible one after the other. Any trailed stars in the 
field of the guided comet are then photometrically Irrelevant. The net effect 
is to have one's cake and eat it, for there will then exist a direct comparison 
between guided photometric standards and the guided comet. The assumptions are 
the same as those that obtain conventionally: that (i) the detector Is uniform, 
and (ii) sky brightness and environmental conditions do not change from one 
part of the procedure to the next. Concerning (i), a lack of detector uniform- 
ity manifests itself as random and systematic error. In the conventional tech- 
nique, random error appears as noise in the calibration, and systematic error 
can be ascertained from subsets of standards across the detector. In the new 
technique, the role of random and systematic errors is the same for the stand- 
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ard star exposure, but a systematic error across exposures can only be inferred 
by extrapolation and inference. Again, the need for controls is apparent, and 
is here accomplished by reversing the sequences of exposures, with rigid adher- 
ence to the maxim 'ceteris paribus'. As to environmental conditions (ii), 
their equality in the new technique is easier to accomplish and to monitor, 
since calibration is actually accomplished at the telescope rather than in the 
laboratory. 

Clearly, these options require a well-designed and executed observing 
program, but this is expected. 


VIII. POLYNOMIAL PROPERTIES 
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§8.1 BAKER-SAMPSON POLYNOMIALS 

The standard rectification function of Baker (1926) and Sampson (1926) 
has the form (10°-1) P , where D is the specular photographic density and p is 
positive. The generalized polynomial expansion is: 

(8.1.1) P = I C n (10 D -l) pn . 

n=l 

At least one coefficient must be positive in order that P be positive in any 
specified domain of interest of (D,logP). Any component P ^ (10 -1) P tends to 
zero as D pn . It is the low numbered terms in the expansion that govern the low 
D asymptote and thus emulate the low density part of the H&D curve. Evidently, 
the lowest numbered coefficient must be positive for real exposures. 

At large D, the n term increases as: 

(8.1.2) P % 10 Dpn , 
so that: 


(8.1.3) dD/dlogP = 1/pn . 

This is the slope of the high D asymptote, which is less for greater p and n. 
Thus for a given p in equation (8.1.1), the last term in the expansion: 

c N (io D -i) pN 


governs the final slope of the polynomial expansion, regardless of the size of 

C . In fact, if C„ is negative, P may eventually turn negative, and its log- 
N N 

arithm is then imaginary. This is the opposite of what is required to model 
saturation and the approach to solarization. 

The condition for a maximum in logP vs. D may be derived by considering 
the competition between only two terms, viz. the last one of the expansion and 
any antecedent term, since it is the last term which will ultimately cause the 
modelled H&D curve to 'turn over' while the antecedent term attempts to model 
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the low density asymptote correctly. Let the terms have expansion numbers m 
and n (m < n) . Physical meaning obtains only when P > 0 and D > 0; i.e.: 
(8.1.4) P = C m (10 D -l) pn + C n (10 D -l) pn > 0, (D > 0). 

A maximum will exist in logP when dlogP/dD = 0; but since P is itself 
never zero or infinite in the domain of interest, the condition is dP/dD = 0. 


Thus: 


(8.1.5) 


mC (10 D -l) pm + nC (10 D -l) pn 
m n 


Dividing equation (8.1.5) by m and subtracting from equation (8.1.4) gives: 


( 8 . 1 . 6 ) 


(1 - -) C (10 D -l) pn > 0 


Thus C < 0 since n > m, and the condition for a real maximum in logP is: 
n 

(8.1.7) p > 0, n > m, D > 0, P > 0; > 0, C n < 0 . 

Condition (8.1.7) alone is insufficient to guarantee a respectable 
range of D for which logP is real; normally one expects to model a range of 
densities up to, say, 3 or 4. This imposes a further restriction on the mag- 


nitude of the coefficients. The value D at which a maximum in logP will 

max 

occur is governed by equation (8.1.5), which can be rewritten as: 


( 8 . 1 . 8 ) 


- 1 + (-mC /nC ) 1/p(n ~" ,) . 


But any respectable value of D must be such that 1 makes a negligible con- 
tribution to the r.h.s. There is therefore a practical constraint on the ratio 
of coefficients given by: 

(8.1.9) p(n-m)D ^ log(-nC /mC ), (D » 1, C < 0, n > m) . 

For example, if m=l, n=2, p=0.25, D = 4, equation (8.1.9) requires that: 

I C 2 | * C^/20. 

In general, in the event that the last coefficient of the expansion is 
negative, and In order for a substantial area of the (D > 0,logP) plane to be 


occupied by real solutions: 


( 8 . 1 . 10 ) 


<< C , (n > m, C < 0). 
m n 
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§8.2 SATURATION AND INFLECTION 

When all coefficients in equation (8.1.1) are positive, then the entire 
(D > 0,logP) plane is occupied by real solutions. The resulting curve is wobbly, 
but the simplest wobble is an inflection that results from the competition be- 
tween two terms of the series, say m and n < m. Results of the previous section 
show that the lower order term governs the low D solution, and the higher order 
term governs the high D solution; the low order coefficient must always be 
positive, and when the higher order coefficient is negative It produces a 
maximum in logP. But when the higher order coefficient is positive, it pro- 
duces a slope of 1/pn which is less than the slope 1/pm. Consequently, inflec- 
tion will always occur as a flattening with increasing D, and this is now 
ultimately the correct form by which to account for saturation. 

The condition for inflection is d 2 logP/dD 2 = 0. For two terms, P is 
given by equation (8.1.4), viz.: 

(8.2.1) P = C_(10 D -l) pm + C„(10 D -l) pn 


m 


n 


Thus points of inflection occur whenever D satisfies: 
[m + nC(10 D -l) p(n_m) ] [1 + C(10 D -l) P(n_m) ; 


( 8 . 2 . 2 ) 


P(n-m) 


c = c /c . 

n m 


= pC (n-m) 2 10" (lCT-l) 1 

where : 

(8.2.3) 

Condition (8.2.2) is better written as a quadratic in C: 

(8.2.4) m + C[m+n-p(n-m) 2 10 D ](10 D -l) P(n ' m) + C 2 n(10 D -l) 2p(n_m) - 0. 

It suffices merely to consider the case m-1 , n=2. Then equation (8.2.4) is: 

(8.2.5) 1 + C[(3-pl0 D )(10 D -l) P ] + C 2 [2(10 D -l) 2p ] = 0, 
which has roots: 

C = [pl0 D -3±A-6pl0 D +p 2 10 2D ]/4(10 D -l) P . 


( 8 . 2 . 6 ) 
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These roots are real whenever: 

(8.2.7) P 2 10 2D - 6pl0° +120. 

Let Y be plO D . Then: 

(8.2.8) [Y-(3+2< / 2)] [Y-( 3-2^2) ] 2 o . 

For the inequality to hold, each factor must have the same sign; either: 

Y-( 3+2^2) 2 o and Y-(3-2» / 2) 2 o, 

or: 

Y-(3+2' / 2) ^ 0 and Y-(3-2' / 2) * 0. 

In other words: 

(8.2.9) 3-2' / 2 2 y 2 3+2> / 2 . 

Consider each limit in turn. 

For the upper limit Y - 3 - 2^2 of equation (8.2.9), it is easy to show 
that C is always negative regardless of the sign of the square root in equation 

(8.2.6) . Thus: 

(8.2.10) plO° ^ 3 - 2^2 , 

corresponds always to negative C. From the previous section, these are path- 
ological H&D curves that first show inflection, before ultimately reaching a 
maximum in logP and reversing course. 

For the lower limit Y - 3+2^2 of equation (8.2.9), regardless of sign 
in equation (8.2.6), C is always positive. Thus: 

(8.2.11) plO D = 3+2^2 , 

is the more desirable case, since inflection bends the solution in direction of 
increasing logP. Two roots for C are ensured by the two signs in equation 

(8.2.6) , except when the equality holds, in which case: 

C = (pld D -3)/4(10 D -l) P . 


( 8 . 2 . 12 ) 
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The array of permissible values of C, and the densities at which they occur, 
are displayed in Table 8-1 for positive C, and in Table 8-2 for negative C. 



Table 8-1 


Values of log C(C=C/C >0, n=m+l=2) for which inflection occurs in a curve 

n in 

modelled by the Baker-Sampson polynomial expansion: C + C (10 D -l) pn . 

ra n 


D 

P = 

0.1 

p = 0.2 

P = 

0.4 

5 

3.1988 

-4.4500 

2.9999 

-5.3009 

2.3010 

-6.6020 

4 

2.2977 

-3.3987 

2.1994 

-4.1004 

1.7007 

-5.2017 

3 

1.3857 

-2.2866 

1.3935 

-2.8944 

1.0979 

-3.7986 

2 

0.3259 

-1.0257 

0.5273 

-1.6265 

0.4683 

-2.3658 

1.7656 

-0.3263 

0.2727 

-1.2770 

0.3015 

-2.0090 

1.5 

- 


-0.1941 

-0.7031 

0.0794 

-1.5693 

1.4645 

- 


-0.4404 

0.0445 

-1.5050 

1.25 

- 


- 

| 

-0.2408 

-1.0401 

1.1635 

- 


- 


-0.6036 

1 

0 

: 


_ 


_ 





Table 8-2 


Values of C (C = C /C < 0, n = ra+1 = 2) for which inflection occurs in a curve 

n m 

modelled by the Baker-Sampson polynomial expansion: C^(10^-l)^ m + C^( 10^-1 . 


D 

P = 

0.1 

P = 

0.15 

P = 

0.17 

0.3 

- 

- 

- 



0.2344 

-0.7312 

- 

- 

- 

- 

0.2 

-0.6777 

-0.8214 

- 

- 

- 

- 

0.1 

-0.6764 

-0.9685 

- 

- 

- 

- 

0.05836 

-0.7023 

-1.0493 

-0.9458 

- 

- 

0.05 

-0.7112 

-1.0707 

-0.9239 

-1.0172 

- 

- 

0.01 

-0.8258 

-1.2843 

-1.1100 

-1.3915 

- 

- 

0.00400 

-0.9036 

-1.4116 

-1.2667 

-1.6084 

-1.5676 

0.001 

-1.0372 

-1.6239 

-1.5556 

-1.9872 

-1.9287 

-2.0435 

0.0001 

-1.3055 

-2.0452 

-2.1958 

-2.8100 

-2.8419 

-3.0353 





IX. NUMERICAL TESTS 
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§9.1 Input Data 

Verification of the theory of Zou, Chen, Peterson (1981) has been 
accomplished by these authors themselves, and by Klinglesmith and Rupp (1984), 
and Speck et al. (1987), using M31. These works were concerned solely with 
derivation of the polynomial coefficients, and did not address the problem of 
the derivation of the extinction coefficients and sky brightness. Walterbos 
and Kennicut (1987) also studied M31, but relied solely on laboratory 
calibration. Warnock and Klinglesmith (1984) have worked on the problem of 
trailed standards. 

As a first step in the development of a data reduction algorithm for 
the general theory of standard star calibrators, it is necessary to formulate a 
problem that will serve to test it and to reveal its numerical properties. We 
seek to answer the general questions: (I) does the theory work In principle, 
and in particular, can the primary extinction coefficient be derived? and (ii) 
are the numerical techniques satisfactory, and what are their properties? 

In order to accomplish this first step, it is necessary to input den- 
sities that are reasonable facsimiles of data to be expected in real cases. 
Since the data available for M31 had standard stars grouped over a very narrow 
range of zenith distance (cf. Speck 1987), and were moreover secured at small 
zenith distance, they were unsuited to a test for a complete solution. With 
the help of sensitometric densities taken from UK Schmidt plates, and from our 
prior analysis of M31, a canonical data-set was established comprising 9 stars 
divided into 2 groups; stars #1 - #4 were located at zenith distance z * 60°, 
and stars #5 - #9 were located at z * 66°. The Input data are summarized 
in Table 9-1. The division into two groups Is necessitated by the fact that 
the primary extinction coefficient cannot be determined unless there is a 
sufficiently large range of z (§6.4). The chosen range of 6° is large enough 
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to eliminate any numerical effects owing to ill-determined matrices, yet 
small enough to correspond to any reasonable range of z accessible to wide- 
field cameras. The background sky + fog contribution was established for each 
group of stars, and the background fog contribution, assumed uniform over the 
field, was similarly fixed. 

Each standard star, sky, and fog frame was comprised of 10 x 10 pixels. 
The number of square arc seconds subtended by a pixel is given by: 

a - (0.001 1 s) a 

where 1 is the pixel length in microns, and s is the plate scale in ”/mm. For 
values of 1 * 20 microns/pixel and s « 100 ”/tnm, a is 4 sq. arc seconds /pixel. 
Star image sizes were determined from the formula of Liller and Liller (1975) 
pertinent to Palomar 1.2m Schmidt plates. In all cases, air masses were simply 
equated to sec z. All input data are therefore well within range of normal 
expectation. 
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§9.2. Tests of the Algorithm. 

Since the primary extinction coefficient is far and away the most 
important one, the primary goal is to see whether there is a reasonable expect- 
ation of its derivation. We describe here the results of tests which are the 
simplest conceivable, yet which are still definitive; viz. tests of the case 
where (i) the tertiary extinction coefficient is ignored, (ii) the secondary 
color coefficient is assumed known, and (iii) only the primary coefficient and 
one or two Baker-Sampson polynomial coefficients are regarded as unknowns. 

Thus the inclusion of the secondary color coefficient and stellar color indices 
(Cl) in Tables 9-2 and 9-3 are for present purposes to be regarded as non- 
participa tory factors insofar as the present tests are concerned, even though 
their values were used to evaluate participatory coefficients [e.g. the l.h.s. 
of equation (7.2.2)]. 

When there is only one polynomial coefficient C, the problem (7.2.2) 
becomes inherently linear in log C and the extinction coefficients, and is thus 
more amenable to direct analytic treatment (§9.4). For the sake of simplicity, 
we describe the parallel numerical implementation of Newton's method as if this 
were the case. 

A program in GWBASIC was written to implement the theory of §7.2. 
Equations (7.2.3a,b) and (7.2.5) were programmed to evaluate the 2x9 matrix A 

r 

of equation (7.2.6), whose 2 columns are necessitated by the 2 unknown correc- 
tions to k^ and the polynomial coefficient C in matrix X, and whose 9 rows 
correspond to the 9 standard stars. Similarly, the 9 rows of matrix Y of 
equation (7.2.6) are evaluated at each iteration for each of the 9 standards. 
Evaluation of the matrix (A A) A proceeds by the usual means (e.g. Natrella 
1963; but note that matrix elements there are opposite to the conventional RC 
- or Row/Column - notation!); this matrix will be merely 2 x 2 if there are 
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only 2 unknown coefficient corrections X, In double-precision, the solution of 
equation (7.2.6) can be found by implementing Cramer's Rule without fear of 
ill-conditioning; our results show that Cramer's Rule works for 2 or 3 un- 
knowns. In the case of more unknowns, more sophisticated methods are necessary 
(e.g. Branham 1988). 

To demonstrate the numerical accuracy and correctness of the program, 
we have evaluated the magnitudes m^ appearing in the terms of matrix Y [i.e. 
the r.h.s. of equation (7.2.5), i = 1,2,. ..,9] such that each term is zero. 

This ensures that the the optimization parameter p, and the coefficients k^ 
and C, used in the calculation, should be the ones that obtain after imple- 
mentation of the iteration algorithm. Double precision Is used throughout, 
with results given in Table 9-2. When the iteration algorithm is then run to 
convergence, the results shown in Figures 9-1 and 9-2, and in Table 9-3, are 
found; the Input data are recovered, so the program works correctly and with 
sufficient accuracy (at least 5 significant figures). Moreover, the test shows 
that the range In z is sufficient to ensure that the determinant A is well- 
conditioned . 

For the test data used here, and for p in the range of about ±0.05 about 
the correct value, Improvement in the coefficients Is found at a rate of an 
order of magnitude or more per Iteration. Farther from the correct value of p, 
convergence Is slower as expected, and Is marked by two requirements: (i) the 
guessed value must be closer to the correct one for convergence to occur, and 
(ii) It is necessary to settle for fewer significant figures before divergence 
sets in. Ultimately, the requirement (i) above becomes absurdly stringent, and 
for all practical purposes the trial solutions diverge. This quantifiable be- 
haviour Is canonical, and is programmable in the general case. 

Owing to the linear properties of the perturbed system, k^ converges 
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at a rate that keeps pace with the polynomial coefficient C, and this means 
that it must be derived formally to an accuracy unwarrented by the observations 
if C is to be accurate to, say, 4 significant figures. It is to be expected 
that a similar condition obtains in the general case, and in fact it might be 
quite satisfactory simply to regard k^ (and perhaps k^) as given. This should 
substantially improve the convergence properties of the algorithm in the 
general case by allowing every bit of information from the standards to 
contribute where their need is greatest, viz. in the determination of k^ and 
the polynomial coefficients. 
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§9.3. Magnitude Errors 

Standard star magnitudes will cutomarily be in error, and these may be 
systematic or random. The effects of random error in a given number of stand- 
ards will be to reduce the ability of the algorithm to converge; this can be 
surmised by carrying the argument to the extreme, when input magnitudes in 
equation (7.2.5) are so poor that the polynomial differences there are in- 
capable of reproducing them. Random error is meliorated by increasing the 
number of standards, but the errors must continue to be randomly distributed in 
zenith distance as well if is to be accurately determined. 

Table 9-3 contains solutions for two worst cases in the case of 9 
standard stars: viz., when the same absolute error of 0.03 mag is systematic- 
ally distributed in z, such that (i) all are too large (i.e. too faint) at 
the smaller z, and too small (i.e. too bright) at the larger z; and (ii) vice 
versa. In the test case here (Tables 9-1, 9-2) the set at smaller z contains 4 
stars, and the set at larger z has 5 stars, so that the median of the random 
errors is as close to zero as needed for this illustration. 

The convergence properties in these two cases are shown in Figure 9-2. 
The upper and lower curves are the envelopes which contain all other solutions 
resulting from any other z-distribution of error (provided the absolute error 
is still 0.03 mag and the choices of z are as stated). In other test cases not 
reported here, convergence is slowed by larger magnitude errors, as expected, 
until finally for large enough errors there is a failure to converge. The num- 
erical results confirm that the algorithm is able to handle situations where, 
owing to observational error or differences in system response, the standard 
star magnitudes and the observed densities are not always precisely compatible. 

The results of Table 9-3 show that when magnitude errors are always 
positive at the lower z and negative at the higher z, the derived extinction 
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coefficient is greater than it would be if there were no errors (and vice 


versa). This may be difficult to see at first glance, but a convincing argu- 
ment proceeds as follows. If magnitudes at the lower z are perceived as pre- 
ferentially too faint, i.e. the magnitude errors are positive, then the ex- 
tinction coefficient would have to be greater to compensate for this. On the 
other hand, the negative magnitude errors which cause the perceived magnitudes 
at the higher z to be too bright, do not require a lower extinction coefficient 
as might be expected, because the air mass is larger there, sufficiently large 
in fact that the negative magnitude error Is accounted for (and vice versa). 
This result is seen also with the help of an analytic treatment. 
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§9.4 Analytic Solutions 

Reduce the system (7.2.2) to its essentials for two unknowns, k and C: 
(9.4.1) m ± + kM ± + 2.5 log C - R ± , (i - 1,2,.. .,1) 

where R^ (i - 1,2,.,,, I) is an abbreviated notation for the polynomial differ- 
ences in equation (7.2.2), and are assumed known. M ■ sec z is the air mass. 
Suppose that equation (9.1.1) is exactly satisfied for all i, i.e. that magni- 
tudes exist in theory which precisely satisfy the system. These "true" magni- 
tudes will presumably differ slightly from the standard star input magnitudes. 
Suppose that these small differences are the same in absolute value for all i, 
but are systematically incurred such that: 

(9.4.2a) m i (new) - m ± + dm, (i - 1,2 Ij, M - M^) 

and 

(9.4.2b) m^(new) - m ± - dm, (i - Ij+1 ,1^+2 , . . . ,1 , M = M 2 ). 

The errors will engender small departures from the "true" coefficients k and C, 
such that : 

(9.4.3) k(new) - k + dk, 
and : 

(9.4.4) logC(new) - logC + dlogC, 

will be the new solution. The r.h.s. of equation (9.4.1) will not change of 
course. Equations (9.4.1) to (9.4.4) give: 

(9.4.5a) m ± + dm + M (k+dk) + 2.5(logC+dlogC) - R^ (i - 1,2,...,^), 

and : 

(9.4.5b) m - dm + H 2 (k+dk) + 2.5(logC+dlogC) - R j ,, (i - 1^1,...,!). 

When there are no magnitude errors, equation (9.4.1) holds, so that: 

(9.4.6a) dm + M^k + 2.5 dlogC - 0, (i - 1,2,...,^), 

and: 

—dm 4* M 2 dk 4“ 2.5 dlogC B 0, (i * X ^41 ,...,!). 


(9.4.6b) 
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Since each system of equations is now independent of i, any one equation from 
each set suffices to determine dk and dlogC. Thus: 

(9.4.7) dk*+ 2dm/(M2“M^), 
and: 

(9.4.8) 2,5 dlogC = - dm^+M^/^-M ) . 

It is easily shown that these expressions account precisely for the results in 
Table 9-3 that have been obtained by numerical implementation of Newton’s 
method , 

Equations (9.4.6a,b) reveal another interesting fact that might not be 
readily apparent. If the signs preceding dm in each equation were to be the 
same, then dk = 0, but dlogC = ±0.4dm. In other words. If there is a truly 
systematic error between the standards and the measured responses, regardless 
of placement in z, then k should still be correctly determined, because k is 
determined differentially. The systematic error is faithfully transmitted to 
the photometric expansion coefficient instead. In practice errors will always 
exist, but if these are normally distributed in z, accurate values for 
extinction should be derivable to first order, provided that there are enough 
standards distributed over a sufficient range of z, and provided that the 
detector is well-behaved for trailed standards. A synopsis of these results 
has been presented elsewhere (Usher 1989). 



Table 9-1 


Summary of Input Densities 
























Table 9-2 


Exact Solution Parameters 
p - kj - 3/10; (k 2 - 0.03); C « 2xi 0~ S 


Star # 

mag 

Cl 

Z 

i 

17.63254.. 

0.5 

60° 

2 

15.93716.. 

1.0 

60° 

3 

15.01066.. 

0.9 

60° 

4 

16.67302.. 

0.7 

60° 

5 

15.71947.. 

1.3 

66° 

6 

14.59827.. 

-0.2 

66° 

7 

17.81131.. 

-0.1 

66° 

8 

16.64795.. 

1.0 

66° 

9 

i 

15.19134.. 

0.8 

66° 




Table 9-3 


Summary of Converged Solutions (k^ ■ 0.03, ■ 0) 


magnitude errors 

P 

k l 

CxlO 8 

W 60O) 


+/-0.03 (z=60®66°) 

0.25 

0.4308 

1.5288 

22.22 


0 (exact) 

i 

3/10 

2 

21.93 



sky 


( 66 °) 


21.93 

21.64 


HU 

fog 

21.13 

20.84 


-/+0.03 (z=60°66°) 


0.25 


0.1692 


20.55 
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X. STANDARD STARS FOR THE P/HALLEY APPARITION 
§10.1 The Search Area. 

The Large-Scale Phenomena Network (LSPN) of the International Halley- 
Watch has provided large numbers of wide-angle images of P/Halley taken with a 
variety of cameras located at sites across the world. For the most part, these 
photographs were taken with instruments that were not equipped with tube spots 
or calibration wedges, so that calibration must be effected through the use of 
field stars of known magnitudes and colors. A list of stars of known identifi- 
cation, position, magnitude, and color, has been compiled for this purpose; the 
entries comprise the Special Halley-watch Identification and Photometric Star 
Catalog (hereafter, the SHIPSCat) and have been given In a separate report to 
the LSPN (Dorband and Usher 1985; hereafter, Report II). 

In selecting the boundaries of the search area. It was necessary to 
account for four mutually dependent effects: (i), the angular extent and 
orientation of the tail; (ii), the sizes of the fields of the network tele- 
scopes; (iii), the Incidence of calibration standards of the requisite accu- 
racy; and (iv), the centers and orientations of the fields of view of the tele- 
scopes along the projected path of the comet. Details of these calculations 
are given in a separate report to the LSPN (Dorband et al. 1985; hereafter. 
Report I). Since it was necessary to compile Report I before the event, it was 
necessary to predict the behaviour of the comet in order to map the avail- 
ability of standard stars along its path. The geometrical aspect and temporal 
evolution of the plasma tail was predicted semi-empirically by Niedner (1984) 
with the help of data from the 1910 apparition. Data on participating tele- 
scopes enabled the concept of "average telescope" to be quantified, by which it 
was surmised that an average f ield-of-view (FOV) might be comfortably taken to 
be 50-100 square degrees. As a result, criteria were established with due 
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consideration for items (i)-(iv) above in order to define the search area. 

These are discussed in Report I. The final search area is shown in Figure 10-1 
and given in greater detail in Report I. 
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§10.2 Properties of the SHIPSCat. 

The question arises as to the incidence of photometric standards that 
are needed to calibrate a plate to the required accuracy. There is a clear 
choice: stars with known photographic magnitudes, or stars with known photo- 
photoelectric magnitudes. The number of photographic standards available per 
unit surface area of the sky can be estimated from the SAO Star Catalog. 
Photographic magnitudes are listed for 131,080 stars. These magnitudes lie 
mainly between 7 and 12. [More than twice as many - almost the full SAO 
catalog - have visual magnitudes.] The average surface density of photo- 
photographic magnitudes is thus just over 3 per square degree, providing 300 
standards over the 10 degree square FOV of the optimal telescope. The accuracy 
of these potential standards is unlikely to be better than 0.1 magnitudes, and 
their usefulness in the implementation of the standard star theory may well be 
problematical . 

On the other hand, the Catalogue of Stellar Identifications (CSI; 
Ochsenbein 1978; Wenger and Ochsenbein 1984) lists about 6300 objects with U, 

B, or V, magnitudes over the Halley path (13% of the sky), for an average 
surface density of about 1.2 per sq. deg. Magnitudes in other systems 
are far fewer; Geneva photometry amounts to 10% of this value, and is statis- 
tically insignificant in the quest for calibration. Stromgren narrow-band mag- 
nitudes are the nearest competitor, but the (b-y), ml, cl, and H-beta indices 
are still only 30% of the broadband magnitudes. Inasmuch as the photographic 
emulsions likely to be used to record the wide-field phenomena are more similar 
to the broadband Johnson filters, we confined the search to U,B, and V. 

The number of B standards available on a 10 degree square plate is thus 
about 100, mainly between magnitudes 5 and 13. Considerable numbers of these 
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may be lost to saturation without hindering a photometric history of P/Halley. 

Difficulties in calibration might be encountered only for FsOV of small 
solid angle , but these are likely to be of less interest to the LSPN. Special 
cases (such as imagery at high z) may require that the number of standards be 
supplemented. This can be accomplished either by generating them ab initio , or 
by generating secondary standards from primary ones. In the first case, high 
photometric accuracy may be achieved relatively quickly on telescopes of modest 
aperture with CCD photometry. In the second case, secondary standards with an 
accuracy of 0.05 mag or better can be generated by the methods outlined. 

A minimum of 10 accurate calibrators In an average FOV are necessary 
to achieve the necessary photometric solution for the large-scale phenomena of 
interest. Barring setbacks unforeseen by the research effort reported here, 
the prospects seem excellent for a quantified photometric and morphological 
history of the latest apparition of P/Halley. 
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